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The  effects  of  subduction  zones  on  local  and  teleseismic  network  data 
for  earthquakes  in  the  Ilorth  Pacific  have  been  deteimined.  Hypocentcr 
location,  fucal  mechanisms  and  slip  vectors  are  all  perturbed  in  ways 
that  can  be  corrected  by  tracing  seismic  rays  through  appropriate  niodels 
of  the  upper  mantle.  The  matching  of  synthetic  seismograms  for  the  first 
30  seconds  of  the  long  period  P signal  with  observations  at  several 
azimuths  yields  depths  of  crustal  events  to  within  10  km  and  relinble 
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estimates  of  focal  mechanisms  and  source  time  functions.  Spectral 
estimates  for  body  waves  based  on  the  autocorrelation  function,  witl 
a Hanning  window,  are  more  reliable  than  those  by  the  FFT 
algorithm.  Station  spectra  can  be  corrected  reasonably  well  for 
source  crust  and  receiver  crust  effects,  but  the  proper  way  to 
model  attenuation  remains  a difficult  problem.  The  use  of  a 
frequency-dependent  modulation  of  T/Q  has  been  tested  and  seems 
promising. 
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Technical  Report  Summary  i 1 

] 

1.  Seismic  ray  tracing  has  been  applied  to  teleselsmlc  and  local 

network  data  from  the  North  Pacific  as  a way  of  examining  the  effects  of  ] 

subductlon  zones  on  seismic  waves.  The*  positions  of  subducting  plates  ] 

can  be  fixed  on  the  basis  of  travel-times  and  amplitudes.  Deeper  than 
normal  hypocenters  fall  In  the  colder*  more  brittle  parts  of  the  plate. 

The  dip  and  azimuth  of  focal  angles  determined  with  plate  effects  taken 
Into  account  are  different  by  large  amounts  from  symmetric  earth  values. 

For  some  distributions  of  stations  with  respect  to  the  source,  slip 

I 

vectors  derived  from  focal  mechanism  solutions  may  be  In  error  by  signi- 
ficant amounts.  Solutions  for  Instantaneous  relative  plate  motions  based 
on  such  data  will  be  adversely  affected.  Only  small  changes  'In  the 
apparent  surface  velocity  and  azimuth  of  the  wavefront  are  predicted 
for  arrays  at  teleselsmlc  distances. 

2.  The  waveform  of  long-period  teleselsmlc  P waves  Is  shaped  by 
numerous  factors.  Including  focal  depth,  focal  mechanism,  the  source 
time  function,  and  the  earth  structure  at  the  source  and  receiver. 

The  sensitivity  to  focal  depth  Is  so  great  that  It  Is  possible  to  fix  j 

] 

the  depth  of  crustal  events  to  at  least  10  km  by  matching  computed  and  < 

observed  waveforms  at  a number  of  stations.  This  resolution  capability  I 

has  been  demonstrated  for  events  in  continental  and  oceanic  crust  in 
the  30*  - 80*  range  of  eplcentral  distances.  ! 

3.  In  addition  to  fixing  focal  depth,  the  technique  of  matching 

Ion ''-period  P waveforms  also  yields  a best— fitting  focal  'inechanlsm  and  a ■ 

source  time  function.  The  technique  has  been  tested  on  12  earthquakes. 

Agreement  with  focal  mechanisms  based  on  standard  techniques  is  very  good  ^ 

in  almost  every  case.  A source  time  function  shaped  as  an  Isoceles 
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triangle  gives  the  best  solutions.  The  method  holds  promise  for  dlscrlm** 
Inatlon  of  explosions  on  the  basis  of  focal  mechanism,  but  the  available 
long-period  P signals  from  moderate  explosions  in  central  Asia  are  Insuf- 
ficient to  permit  a thorough  test. 

A.  Further  studleb  ul  the  methodology  of  computing  source  spectra 
from  digital  seismic  data  have  led  to  the  development  of  a program  based 
on  the  autocorrelation  function  that  produces  more  reliable  estimates  of 
the  spectra  than  those  yielded  by  the  Ft  at  Fourier  Transform.  Experiments 
with  various  data  windows  have  been  conq>l%ted  and  the  Hanning  window  has 
!)Gcn  demonstrated  to  be  the  best.  The  way  to  treat  Inelastic  attenuation, 
or  the  proper  choice  of  a Q-model,  remains  the  least  understood  of  the 
various  corrections  to  observed  spectra  required  by  the  properties  of  the 
real  Earth.  A modification  of  the  standard  practice  of  taking  the  ratio 
of  travel-tlme-to-Q  (T/Q)  to  be  a constant  Is  achieved  by  applying  a 
relaxation  mechanism  for  attenuation,  proposed  by  Solomon.  A frequency  * 
dependent  modulation  of  T/Q,  with  a characteristic  relaxation  time  as 
a parameter,  has  been  used. 

5.  A paper,  "Rupture  Velocity  and  Indication  Efflcienty",  reporting 
work  partly  supported  under  this  grant,  has  been  submitted  to  the  Bulletin 
of  the  Selsmologlcal  Society  of  America.  The  abstract  is  Included  for  part 
of  this  report. 
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1.  Plate  Effects  in  North  Pacific  Subductlon  Zones 


E.  R.  Engdahl  and  Ming-Te  I, in 


Introduction 


Since  Davies  and  McKenzie  (1969)  first  documented  advances  In  P 


residuals,  the  existence  of  hlgh-veluclty  plates  In  North  Pacific  sub' 


duction  zones  has  been 'an?>ly  demonstrated.  In  view  of  the  unrealistic 


nss'imptions  of  a spherically  symmetric  Earth,  the  existence  of  a subduct  ed 


plate  beneath  an.  Island  arc  often  malces  It  difficult  for  geophysicists 


to  study  seismic  sources  in  these  regions. 


Research  Is  now  centering  on  an  improved  understanding  of  upper 


mantle  structure  and  wave  propagation  effects  in  subductlon  zones.  The 


rf  Ccct  s ot  (ironounced  lateral  heterogeneity  in  physical  proper ti>*s  on 


the  seismic  signal  have  been,  discussed  by  numerous  authors  (Jacob,  1972; 


Davies  and  Julian,  1972;  Sleep,  1973;  ToksBz  ct  al..,  1971).  The  spatial. 


distribution  of  intermediate  depth  earthquakes  within  slabs  of  lithosphere 


which  descend  into  the  mantle  is  important  to  our  understanding  of  the 


subductlon  process.  Engdahl  (1973,  1.974)  considered  plate  effects  on 


earthquake,  hypocenters  by  allowing  for  plate  structure  in  the  location 


process  rather  than  assixralng  that  the.  Earth  Is  laterally  homogeneous. 


Because,  of  refraction  effects  the  point  on  the  focal  sphere  correspondin* 


to  the  ray  to  a neismograph  station  will  be  plotted  in  the  wrong  place 


if  that  ray  has  gone  through  a high  velocity  subducted  slab,  but  a 


spherical  Earth  me  del  is  used  to  determine  the  ray  path.  The  size  of  the 


effect  may  depend  on  the  source-slab-station  syst^.  . 


In  this  paper  we  are  concerned  with  effects  of  crust  and  upper 


mantle  structure  on  all  aspects  of  the  interpretation  of  seismic  data 


(travel  tines  and  slownesses,  hypocenter  locations,  wave  amplitudes. 


and  focal  mechanisms).  We  begin  \rlth  thermal  models  and  construct 


velocity  models  of  North  Pacific  subductio!i  /ones.  Wc  apply  seismic  ray 


N.U.  Sleep  is  one  of  the  co-auCliors.  This  paper  has  been  submitted  to 
Bulletin  of  the  Scismological  Society  of  America. 
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tracing  to  study  the  effects  of  these  raodels.  To  provide  constraint  on 


> 


the  analysis  of  focal  parameters  and  focal  mechanism  solutions  we  only 
use  data  from  regions  having  good  local  seismograph  networks.  Finally, 
we  demonstrate  how  these  new  data  enable.us  to  better  interpret  the  sub- 
duct Ion  process. 


F 

■ 
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Velocity  Models 

Niuserlcal  theznal  models  of  slabs  (Toksbz  et  al. , 1973)  were  calculated 
for  kinematic  parameters  relevant  to  the  Kuril,  Aleutian  and  Alaska 
subdue tlon  zones.  Velocity  models  were  constructed  from  the  thermal 

I • 

models  using  a linear  relation  between  velocity  and  temperature  perturbation 
(eqn.  1;  Sleep,  1973)  and  the  Herrin  et  al.  (1968)  velocities  as  the 
unperturbed  base  model.  A thermal  coefficient  of  seismic  velocity  of 

_A 

9 X 10  kos/sec*C  seemed  to  give  velocity  contrasts  in  each  region 
consistent  with  published  values  for  the  Aleutians  (Jacob,  1972;  Sleep, 

1973)  and  for  the  Tonga-Fljl  region  (Mitronovas  and  Isacks,  1971).  We 
purposely  us€id  a base  model  without  a low  velocity  zone  or  gradients  due 
to  phase  changes  to  simplify  the  ray  tracing  and  because  these  complications 
were  not  warranted  by  the  problem.  Tiic  resulting  velocity  models  shown 
in  subsequent  figures  are  characterized  by  high  velocities  in  the  iipper 
half  of  the  plate,  on  the  average  8 - 10%  higher  than  the  surrounding 
mantle,  corresponding  to  the  cooler  temperatures  of  the  thermal  models. 

The  location  of  each  plate  is  given  by  the  distance  of  some  point 
on  the  model  to  the  center  of  curvature  of  the  suhduction  zone.  This 
allows  us  to  project  a I v;o-d imension.il  model  into  three  dl)neri;:ions  by 
assuming  symmetry  along  the  arc  about  the  center  of  curvature.  The 
radius  and  center  of  curvatu-.e  are  determined  by  a least  squares  fit  to 
the  trench  axis  or  to  the  volcanic  arc.  Parameters  fitted  to  the  trench 
axis  of  the  Kurils  and  of  the  central  Aleutians  and  to  the  volcanic  arc 
in  Alaska  are  given  in  Table  I.  All  sections  shown  in  this  paper  are 
projections  with  respect  to  the  distance  to  the  center  of  curvature. 

The  exact  location  of  the  reference  point  on  l e plate  model  is  not 
uniquely  determined,  however,  but  has  to  be  estimated  from  ob.served  travel 
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times  and  amplitude;,  (Jacob,  1972;  Sleep.  1973).  A good  first  guess  Is 
that  the  Benioff  zone  defined  by  conventional  hypocenters  falls  in  the 
higher  velocity  colder  portion  of  the  plate  at  a depth  of  100-125  km 
beneath  the  volcanic  arc  (Engdahl,  1973). 
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Seismic  Ray  Tracing 

The  problem  of  ray  tracing  in  a generally  heterogeneous  medium  has 
been  treated  using  the  calculus  of  variations  and  Fermat's  principle 
of  stationary  time  (Julian,  1970).  The  solution  In  geocentric  spherical 
polar  coordinates  Is  a system  of  first-order  differential  equations  giving 
the  variation  with  time  of  the  position  (r,  0,  Q)  of  a point  on  a ray 
and  the  direction  of  propagation  (1  ° Incident  angle,  o ^ azimuth)  from 
the  point  in  terms  of  the  wave  speed  (v)  and  Its  spatial  derivatives  in 
the  medium.  Julian  has  developed  a computer  program  to  treat  this  gen- 
eralized formulation  and  this  program  forms  the  basis  for  the  calculations 
prt.'.S',aLed  in  this  paper.  Only  the  size  of  computer  core  limits  the  delnil 
to  which  V can  be  represented  In  the  model.  To  study  subductlon  zones 
it.  is  ronvenient  to  represent  the  structure  by  a two-dimensional  grid, 
uonual  to  the  arc,  in  which  velocities  are  specified  at  discrete  points. 

With  this  representation  any  degree  of  structural  complexity  may  be 
defined  subject  only  to  the  grid  interval.  The  wave  speed  and  Its  derivatives 
may  then  be  determined  numerically  at  any  point  within  the  structure 
using  cubic  spline  Interpolation.  In  practice,  time  steps  are  chosen 
so  tlut,  as  the  ray  propagates,  the  .interpolation  grids  along  the  ray 
overlap.  The  grid  size  is  also  chosen  to  be  about  the  same  order  as  the 
wavelength  of  observed  data.  The  velocity  models  shown  in  Figures  2, 

5 and  9 are  represented  by  grids  at  10  km  depth  Intervals  and  10,  12.8 
and  15.07  km  surface  distance  Intervals  to  the  centers  of  curvature, 
respectively. 

The  next  step  is  to  incoiiwrate  the  plate  model  in  the  hypocenter 
det*?.rmi nation,  using  seismic  ray  tracing  in  the  manner  described  by 
Engdahl  (1973).  This  method  differs  from  the  conventional  formulation 
only  in  that  the  coefficients  of  the  station  condition  equatinns  are  now 
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explicitly  knotm  from  the  ray  vectors  to  each  station  from  the.  source. 


The  problem  of  finding  that  particular  ray  that  arrives  at  a given  station 

is  formidable,  however,  and  is  now  done  by  iteratively  shooting  rays 

« 

until  we  find  one  closer  than  1 kn  to  the  desired  end  point.  The  wavefront 
of  tliat  ray  is  then  propagated  the  small  distance  remaining. 


9 

Data 

To  reduce  computational  expense  the  events  studied  In  each  subductlon 
zone  were  relocated  using  local  stations  and  a fixed  teleselsirlc  network 
of  UWSSN  and  selected  Russian  and  CanadlaV  stations.  P onsets,  first 
notions  and  amplitudes  were  read  from  70  nm  film  copies  of  the  records  of 
WWSSN  stations.  Data  for  local  stations  in  the  Kurils  were  taken  from 

, t 

bulletins;  data  from  the  Amchitka  network  (Engdahl,  1973)  were  read  from 
Develocorder  film  and  tape  playbacks;  and  data  from  the  USGS  Alarkan 
network  were  kindly  provided  by  John  Lahr. 

Station  corrections  developed  by  Sengupta  and  Julian  (1976)  from  a 
worldwide  distribution  of  deep  focus  events  were  applied  to  those  arrivals 
beyond  30*  distance.  '®'ocal  mechanisms  were  determined  from  long  period 
P-wavc  first  motions  on  TWSSN  records  whenever  possible  and  from  short- 
period  first  motions  at  other  key  stations. 


I 
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Kuril  Islands  ' 

Eleven  earthquakes,  large,  enough  to  be  well  recorded  on  long-period 

WWSSN  instruments,  were  chosen  for  special  study  In  the  Kuril  Islands 
(Table  II).  These  events  ave  distributed  along  the  arc  and  with  depth 
(Figures  1 and  2)  to  insure  a good  sampling  of  the  structure.  Several 
closely-located  very  deep  events  were  used  In  the  hope  of  determining 
station  corrections  apprOijriate  to  the  Kurils.  This  was  not  possible 
with  so  few  events  because  the  residuals  did  not  appear  to  be  slowly 
varying  functions  of  the  spatial  coordinates.  This  phenomena  has  also  been 
reported  by  Engdahl  (1975)  for  deep  focus  Fiji— Tonga  events. 

Nine  of  the  events  were  relocated.  In  the  manner  previously  described, 
using  seismic  ray  tracing,  and  a long  slab  model  of  the  Kurils  (Figure  2) . 
The  position  of  the  model  was  fixed  by  locating  the  volcanic  arc  .'hU  C. 

The  two  shallow  events  in  the  left  hand  comer  of  Figure  2 were  not 
relocated  but  were  selected  for  a special  study  of  plate  effects  on  slip 
vector  determinations  to  be  discussed  in  a later  section. 

The  model  is  undoubtably  a gross  oversimplification  of  the  overall 
arc  structure  but  is  useful  in  examining  plate  effects  on  P-wave  residuals, 
slownesses,  focal  mechanism  determinations  and  long-period  amplitudes. 

In  addition,  seismic  ray  tracing  allows  local,  regional  and  teleselsmlc 
data  to  be  used  jointly  without  the  disadvantage  of  huge  biases  produced 
by  ray  paths  travelling  primarily  in  the  high  velocity  plate,  especially 
to  local  stations  (Mltronovas  and  Isacks,  1971). 

For  the  deepest  events  in  this  set,  travel  time  advances  of  up  to 
7 sec  to  stations  along  the  arc  are  observed.  Travel  times  to  teleselsmlc 
stations  from  intermediate  depth  events  were  commonly  decreased  by  up  to 
3 sec.  The  joint  effect  on  location,  compared  to  symmetric  earth 
calculations  is  to  shift  the  hypocenters  to  the  northwest  and  deeper. 


11 


Only  small  changes  of  up  to  0.05  secAltg  in  slowness  and  0.5“  in  azimuth 
of  the  wavefront  are  predicted  for  arrays  at  teleseisraic  distances 
(>  30“). 

Plate  structure  may  also  severely  effect  the  apparent  location  of  a 

« • 

ray  .lonvlnf,  the  focus  and,  hence,  in  some  instances  produce  significant 
changes  in  l:ocal  meclianism  determinations.  Focal  angles  for  the  nine 
relocated  events  were  found  to  differ  from  synmetric  earth  values  by  up  to 
20*  in  take-off  angle  and  35“  in  azimuth  for  teleseijjms.  Focal  mechanisms 
were  determined  for  all  the  relocated  hypocenters  using  only  traced  rays. 

All  but  two  of  these  events  were  also  f.tudied  by  Stauder  and  Mualchin 
ii.‘.  Lug  long-period  P waves  and  S wave  polarization  angles  with  a 
symmetric  earth  model . The  traced  solutions  have  nodal  planes  that 
generally  conform  ft*  those  of  Stauder  and  Mualchin  but  with  small  de- 
I'actutc;'  :.ti  dip  ami  strike  of  up  to  and  10“,  respectively. 

Although  w'  have  calculated  the.  details  of  the  velocity  distribution 
from  i henna  I contours,  the  plate  is  still  not  uniquely  located  in  the 
Earth,  Sleep  (197  3)  has  shown  that  in  the  case  of  the  LONRSHOT  explosion 
in  the  Central  Aleutians  the  amplitudes  arc  particularly  sensitive  to  plate 
position.  Wo  have  chosen  position  C in  Figure.  7.  for  our  starting  model, 
since  this  gives  a depth  to  the  Benioff  zone  beneath  the  volcanic  arc  of 
roughly  125  km  as  suggested  by  Engdahl  (1973).  We  use  essentially  the 
same  technique  as  Sleep  (1973)  to  examine  the  amplitude  of  the  first 
peak-to-peak  motion  recorded  on  long-period  WWSSN  instruments  from  nine 
of  the  Kuril  events.  This  technique  requires  correcting  observed  amplitudes 
for  the  source  mechanism  determined  from  traced  rays  and  for  epicentral 
dist.uicc.  Source  corrections  for  rays  near  nodal  planes  were  icslrictcd 
to  a maximum  factor  of  20.  A combLued  station  and  distance  correction  was 
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obfained  by  using  the  deepest  event  on  30  Aug  1970  to  correct  the 
relative  amplitudes  of  o^’hcr  events.  In  this  way  we  hoped  to  islolatc 
1 he  purely  plate  effects.  Peilods  of  first  arrivals  were  also  measured 
and  the  a.Tiplltudes  corrected  for  Instrument  response. 

Station  corrected  relative  long-period  amplitudes  given  on  a log^ 
scale  for  a Kuril  event  at  a depth  of  114  loa  are  plotted  on  an  equidistant 
azimuthal  projection  about  the  focus  in  Figure  3.  To  demonstrate  the 
distortion  of  ray  paths  for  plfte  positions  A-D  , rays  corresponding  to 
10*  distance  and  30*  azimuth  increments  on  a synmetrlc, earth  were. traced  through 
the  plate  model  and  also  plotted  in  the  same  projection  about  the  focus. 

For  a symraetric  earth  or  a very  deep  earthquake,  this  ray  density  pattern, 
would  be  cqui-spaccd  concentric  circles  and  radii.  This  plot  is  useful 
in  a qualitative  w.iy  since  ciot'.cly  spaced  lines  .suggest  higher  amplitudes 
while  loige  open  .'iveas  suggest  shadow  zones. 

A large  shadow  zone  to  the  northwest  is  predicted  for  this  event  for 
all  plate  positions.  An  analytical  solution  for  wave  propagntinn  in  a 
one-dimensional  slab  predicts  that  long-period  seismograms  In  the  ray- 
theoretical  shadovA  zone  will  typically  have  amplitudes  507.  (jnx  greater) 
of  tho  direct  P-A/ave  amplitude  and  exhibit  wave,  form  broadening  (Ward  and 
Aki,  1975).  The  rcisiilts  plotted  in  Figure  3 suggest  .a  long-period 
.implitude  reduction  by  at  least  a factor  of  2.  There  also  appeared  to  he. 
no  consistent  pattern  of  wave  form  broadening  In  the  period  measurements, 
due  in  part  to  the  fact  that  there  are  unknown  frequency  dependent  effects 
beneath  stations.  A disappointing  result  in  this  study  is  that  for  all 
these  events  the  long-period  amplitude  ray  density  pattern  is  not  sensitive 
enough  to  plate,  position  for  intermediate  depth  events. 
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Central  Aleutians 

Thirty-nine  intermediate  depth  earthquakes  of  magnitude  A or  greater 
were  selected  for  special  study  in  the  Central  Aleutians.  The  distribution 
of  these  earthquakes  and  the  local  seismograph  network  arc  shown  in 

I • 

Figure  A.  Events  less  than  about  125  km  in  depth  are  uniformly  distrihutecl 
along  tiie  arc.  Deeper  earthquakes  are  clustered  behind  an  active  volcano 
(flemisopochnoi)  along  an  extension  of  fiomrs  Ridge,  a large  arcuate 
aseismic  structure  extending  into  the  Bering  Sea. 

The  seven  events  Indicated  by  the  solid  circles  In  Figure  4 and  listed 
in  Table  IT  were  also  detected  teleseismlcally.  In  Figure  5 is  shown  a 
vcloc  ».t.y  model  of  i.he  Aleutian,  slab  oonsLcuctcd  from.  tUenn.'il  contours. 

The  position  of  t.U.s  model  'is  well  constrained  by  numerous  seismic  travel- 
tiiif*  residuals  frc.T  the  I.ONGSHOT  explosion  (Jacob,  1972;  Sleep,  1973)  and 
'•y  astorlnted  amplitude  anomalies  (SLotp,  1973).  This  slab  model  was 
used  to  relocate  .•^Jx  of  the  telcscismic  events,  shown  in  section  in 
Figure  5,  using  .seismic  ray  tracing.  The  shallow  event  was  not  relocated, 
but  is  the  location  of  the  1965  Rat  Islands  mainshock.  The  result 
confirms  an  earlier  conclusion  by  Engdahl  (1973),  using  a more  simplified 
slab  model,  that  the  Benioff  ^one  at  depth  is  comparatively  thin  and 
corresponds  closely  to  the  colder,  more  brittle  region  of  the  slab 
suggested  by  thermal  models. 

Figure  6 is  a section  showing  39  events  that  have  been  relocated 
using  the  same  model  and  seismic  ray  tracing  but  only  with  data  from  the 
local  network.  The  solid  circles  are  the  new  locations  of  the  same 
teleseismic  events  shown  in  Figure  5.  The  relocated  hypocenters  define 
a thin  zone  in  the  upper  half  of  the  plate  corresponding  to  the  region  of 
largest  velocity  (or  thermal)  gradient.  The  thinness  of  this  seismic 


zone  we  believe  to  be  model  Independent.  The  relative  location  within 
the  model,  however,  is  probably  not  well  determined  since  little  Is  laiown 
about  the  details  of  upper  mantle  structure  in  the  wedge  traversed  by 
rays  to  the  local  network. 

With  these  the  best  estimates  of  our  local  Aleutian  events  we  can 
now  hope  to  obtain  more  consistent  composite  focal  mechanisms  deterrain- 
Htions.  If  Intermediate  depth  earthquakes  occur  in  the  col.der  more  brittle 
region  of  the  slab,  we  would  expect  focal  mechanism  solutlon«>  to  suggest 
that  the  descending  lithosphere  acts  as  a stress  guide  that  aligns  the 
earthquake  generating  stresses  along  the  dip  of  the  slab  (Isacks  and 
Molnor,  1971).  I'igute  7 shows  two  classes  of  mechanisms  obtained  in  the 
Central  Aleutians.  Intermediate  events  from  about  70  to  127  km  in  depth 
are  inlformly  characterized  by  dovm-dip  pressure  axes  or  lateral  coctenslon 
along  the  strike  of  the  arc.  Among  the  deeper  events  there  is  a large 
subset  that  appear  to  have  near  vertical  tension  axes.  The  sense  and 
location  of  this  change  in  focal  mechanism  is  very  similar  to  the  results 
of  Stauder  and  Mualch.in  (1975)  in  the  Kuril  Islands^  They  found  that 
intermediate  depth  foci  separate  into  those  v>hich  represent  an  axial 
comjircssion  with  respect  to  the  plate  at  depths  less  than  120-3  AO  km 
and  those  which  represent  an  axial  tension  at  greater  depths.  Three 
previously  published  solutions  for  deeper  than  normal  events  in  the 
Aleutians  (Stauder,  1968;  Isacks  and  Holnar,  1971)  led  Isacks  and  Molnar 
to  conclude  that  axial  tension  was  the  operative  stress  down  the  Aleutian 
slab.  These  events  are  exceptions  to  the  results  of  this  paper  which 
suggest  that  the  slab  is  under  compression  not  tension.  This  would  occur 
if  the  subduction  rate  is  faster  than  the  rate  at  which  the  lithosphere 
is  sinking  (Isacks  and  Molnar,  1971;  Smith  and  ToksOz,  1972). 


Alaskt\ 

The  uses  southern  Alaska  seismograph  data  provide  excellent  data 
for  resolving  the  configuration  of  the  Benloff  zone  In  the  vicinity  of  - 
Cook  Inlet  (Figure  8).  Lahr  et  al.  (1974)  find. an  apparent  thickness  of 
tlie  dipping  seismic  zone  of  only  about  15  kn  when  projected  for  300  km 
along  strlka.  Seven  large  events,  well  recorded  on  the  local  network  rnd 
teleselsmlcally , were  selected  for  relocation  (Table  II) . We  located 
the  plate  by  placing  the  Benloff  zone  at  a depth  of  115  km  beneath  the 
volcanic  arc.  The  plate  model  differs  from  those  previously  studied 
In  that  the  seaward  extension  is  characterized  by  a long  gentle  slope 
(Figure  9) . The  relocated  hypocenters  and  the  1964  Alaska  mainshock 
are  shown  in  section  with  the  model  In  Figure  9.  At  least  two  events 
appear  to  be  located  at  the  top  of  the  slab.  It  is  difficult  to  interpret 
these  results  since  the  relocated  hypocenters  may  reflect  the  effect  of 
upper  mantle  structure  on  the  preponderance  of  local  arrivals  and  the 
possible  rolslocatlon  of  the  slab.  Further  study  of  these  data  wil],  be 


necessary. 
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Slip  Vectors 

It  has  been  shovm  in  this  paper  and  elsewhere  (Toksbz  et  al..  1971) 
that  the  structure  of  subducted  plates  bends  rays  leaving  island-arc 
earthquakes  and  distorts  the  fault-plane  solution  Inferred  from  first 
motions.  Hence,  under  certain  circumstances  of  station  geometry  with  respect 
to  the  seismic  source,  the  strike  of  either  one  or  both  of  the  nodal 
planes  determined  with  s synnetrlc  Earth  model  may  be  significantly 
in  error  in  the  presence  of  a subducted  plate.  In  some  arcs  a systematic 
error  in  the  slip  vectors  may  be  Introduced,  producing  some  difficulty 
In  modelling  relative  plate  motions  (Minster  et  al. , 1977»).  We  have 
chosen  five  North  Pacific  events,  corresponding  to  overthrusting  of  the 
Island  block  relative  to  the  oceanic  block,  for  a special  study  of  this 
effect.  This  Includes  two  Kuril  events  (Stauder  and  Bollinger,  1966b), 
the  Rat  Islands  malnstock  (Stauder,  1968a)  and  the  1964  Alaska  earthquake 
(Stauder  and  Bollinger,  1966a).  The  focal  parameters  for  these  events 
are  listed  in  Table  II  and  their  local  ion  in  section  shown  in  Figures 
2,  5 and  9. 

The  original  P-wave  first  motion  data  and  nodal  planes  published  by 
Stauder  for  these  thrust  earthquakes  are  plotted,  using  the  Herrin  ot  al.. 
(1968)  symmetric  Earth  model,  on  equal-area  azimuthal  projections  of  the 
lower  hemisphere  in.  the  left  hand  column  of  Figures  10  and  11.  The  same 
rays  were  traced  through  the  appropriate  plate  structures  and  the  new  ray 
projections  plotted  with  the  original  nodal  planes  in  the  second  column 
of  these  figures.  The  third  column  shows  the  ray  distortion  produced  by 
the  plate  for  each  event  in  the  manner  previously  described. 

For  the  two  Kurile  events  there  is  characteristically  one  nodal 
plane  dipping  steeply  to  the  southeast  which  seems  to  be  well  determined  by 
the  Cirst  motion  data  and  is  little  effected  by  plate  structure.  The 
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second  nodal  plane  Is  not  so  %rell  defined  but-  Is  constrained  to  dip  f'ently 
to  the  northwest  by  S wave  polarization  angles*  Stauder  and  Maulchln 
(1975)  found  in  many  Instances  that  the  S polarization  lines  do  nri  converf.e 
well  to  the  T axis  but  seem  to  Indicate  an  axis  with  a plunge  some  5 to 
■JO*  steeper  than  the  one  which  bisects  the  dihedral  angle  between  the  two 
nodal  planes.  This  Is  exactly  the  kind  of  deflection  we  see  in  Figure  10 
for  rays  to  stations  northwest  of  the  source.  Since  the  T axis  of  their 
solutions  is  still  consistent  with  the  nodal  planes  we  conclude  that  their 
published  estliaates  wlU  not  be  greatly  changed  by  plate  structure  and 
that  faulting  characteristic  of  normal  depth  earthqiiakes  in  the  Kurils  is 
consistent  with  the  stated  convergence  of  the  oceanic  plate  In  the  direction 

N55*W  with  respect  to  the  continental  plate. 

For  the  Rat  Islands  nainshock  shown  at  the  top  of  Figure  .11  introduction 
of  plate  structure  has  produced  changes  in  the  focal  angles  of  up  to  8 
in  dip  and  22*  in  azimuth  to  some  stations.  The  largest  changes  occur 
for  rays  departing  nearly  along  the  strike  of  the  arc.  Ray  projections  to 
North  America  appear  to  have  undergone  a net  outward  shift  and  counter- 
clockwisa  rotation,  as  these  rays  seek  the  higher  velocity  path  through  the 
plate.  Ordinarily,  these  changes  might  not  affect  the  location  of  the  nodal 
planes,  as  in  the  case  of  Intermediate  depth  earthquakes  in  the  central 
Aleutians  (Engdahl,  1974).  In  this  particular  case,  however,  the  strike 
of  the  auxiliary  plane  is  highly  dependent  on  stations  to  North  America 
(NE  quadrant).  An  auxiliary  plane  with  a 15*  change  in  dip  direction  would 
fit  t\m  data  equally  well  and  illustrates  the  potential  danger.  The  slip 
vector  could  undergo  a corresponding  rotation  from  N39*W  to  N54*W.  The 
1964  Alaska  earthquake  shown  in  the  bottom  of  this  figure  occurs  too  far 
south  of  the  deep  plate  structure  in  Alaska  for  rays  to  sense  the  velocity 
anomaly,  so  that  the  focal  mechanism  rcroalus  well  determined. 
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These  results  have  Important  consequences  for  the  modelling  of 
relative  plate  notions.  Because  of  the  similarity  In  source-station  geo- 
metry along  the  vhole  of  the  Aleutian  arc,  we  could  expect  to  have  the 
same  discrepancy  In  other  solutions  published  by  Stauder.  Minster  ct  al. 
(1974)  relied  entirely  on  North  Pacific  data  to  propose  the  existence 
of  a Bering  Plate,  moving  with  respect  to  North  America.  Traditionally, 
the  lithosphere  north  of  the  Aleutian  arc  has  been  considered  part  of  the 
North  American  plate  (McKenzie  and  Paricer,  1967). 

The  Bering  plate  was  introduced  by  Minster  et  al.  because  of  a 
systematic  misfit  to  the  slip  vectors  of  earthquakes  along  the  Aleutian 
and  Kurile  arc  for  a pole  (NOAM-PCFC)  determined  by  assuming  the  lithosphere 
north  of  the  Aleutian  arc  as  part  of  the  North  American  plate.  The  Alaska 
earthquake,  the  least  effected  by  slab  structure,  does  not  agree  with 


the  BERI-PCFC  pole  and  would  agree  better  with  the  NOAM-PCFC  pole.  The 
Aleutian  earthquakes  had  misfits  to  the  slip  vectors  of  the  same  magnitude 
and  direction  of  the  changes  in  the.  slip  vector  showr  for  the  Rat  Inlands 
earthquake.  The  two  Kuril  earthquakes  used  were  less  than  the  char- 
acteristic slip  vector  for  the  Kprlles  (N55“W)  and  In  the  direction  of 


the  misfit.  Since  the  data  do  not  overwhelmingly  support  the  existence 
of  a Bering  plot,  we  must  conclude  that  there  is  still  considerable 
uncertainty  in  the  need  lor  a Bering  plate  in  the  modelling  of  relative 


plate  motions  in  the  North  Pacific. 
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Discussion 

The  relocated  hypocenters  plotted  in  preceding  figures  are  partly 
model  dependent,  but  represent  nt  least  a first  order  Improvement  in 
absolute  location  over  previous  studies  In  which  plate  effects  have  not 
been  accounted  for.  In  the  Aleutians,  for  example,  source  bias  due  to 
plate  structure  may  be  expected  to  produce  serious  errors,  of  up  to  0.5 
degree  in  teleseismlc  epicenter  determinations  (Herrin  and  Taggart, 

1968).  In  Alaska,  Lahr  et  al.  (1974)  found  using  only  local  data  systematic 
shifts  In  epicenters  of  up  to  25  km  in  distance  and  10  Ion  In  depth,  the 
location  bias  Increasing  with  depth.  Clearly,  if  we  are  to  use  seismicity 
to  study  mechanical  processes  within  subducted  plates,  we  must  have  a 
better  understanding  of  where  the  earthquakes  are  actually  occurring. 

Conventional  hypocenter  locations  in  the  Tonga  Island  arc  (Rkyc.s  et.  al., 

1969;  Mitronovaa  ot  al.,  1969)  suggest  an  upper  limit  to  the  Lhicljiess 
of  the  seismic  zone  of  about  25  km,  but  considering  errors  of  location 
it  could  be  even  thinner.  Wyss  (1973)  reports  from  the  same  dimensions 
of  large  intermediate  depth  earthquakes  In  the  Tonga  region,  that  the 
thickness  of  the  seismic  zone  is  more  nearly  11  km,  Lahr  et  al.  (197A) 
used  data  from  the  southern  Alaska  seismograph  network  (Figure  8)  to  \ 

estimate  an  apparent  thickness  of  the  Benioff  zone  of  about  15  km  when 
projected  for  300  km  along  strike.  Subdivision  of  the  zone  into  segments 
extending  50  km  along  strike  Indicates  that  in  many  areas  earthquakes  ora 
concentrated  within  a zone  less  than  5 km  thick.  Billlngton  and  Isacks 
(197.5)  used  the  method  of  Joint  Hypocenter  Determination  to  relocate 
500  to  600  depth  earthquakes  in  the  northern  part  of  the  Tonga  arc.  They 
found  evidence  for  fault  planes  within  the  relatively  thick  lithospheric 
slab.  This  study  lends  further  support  to  this  evidence  for  Interior 
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deformation  in  subducted  plates.  A remaining  uncertainty  is  the  question 
of  lateral  variation  of  slab  properties  or  for  tearing  at  depth  (Isacks 
and  Molnar,  1971;  Abe,  1972).  The  answer  to  this  question  is  being 
sought  in  a very  detailed  study  of  Central  Aleutian  earthquakes  now  in 

progress. 

Intemedlat.  depth  foci  In  the  KurU.  ere  charecterlted  by  two 
properties  (Steuder  end  Maulchln,  1976):  1)  the  principal  axis  syatei. 

corresponding  to  the  focal  mechanls..  .11  have  the  B axis  nearly  horizontal 
and  lying  parallej.  to  (l.e.  wlthljJ  the  henioft  rone,  and  2)  one  of  the 


other  two  axes,  either  the  P or  T axis  la  aligned  doim  dip  to  the 
descending  plate.  For  loci  above  120  In.  the  earthquakes  are  down  dip 
conpresslon  and  for  foci  deeper  than  140  to  down  dip  tension.  From 
120-140  to  In  depth  both  types  of  foci  are  found  in  close  spatial  re- 
lationship. Velth  C1974)  has  suggested  that  this  sharp  zoning  may  he 
associated  with  the  depth  at  which  the  ollvlng-splnel  phase  transition 
occurs  within  the  Interior  of  the  descending  plate.  He  nuggests  that  the 
volume  change  assocUted  with  phase  transition  could  result  In  down- 
dip  tenslonal  stress  reaching  a critical  value  within  the  Interior  of  the 
Buhducted  plate  :;hUo  down-dip  compreaslonal  atross  could  equally  reach 
critical  values  across  the  phase  transition  boundary,  near  the  outer 
surface  of  the  plate.  The  region  of  ollvlne-splnel  phase  transition  is 
usually  estimated  to  be  250  to  (Turcotte  and  Schubert,  1971)  to  300  to 
(Rlngwood,  1972).  One  or  more  other  phase  changes,  such  as  basalt  to 
eclogue  and  pyrolitc  to  garnet  pyrollte  are  more  likely  to  take  place 
in  the  upper  200  to  of  the  mantle  and  these  phase  houndarles  may  he  marked 
by  complications  ..f  the  low  velocity  zone.  Those  phase  changes  were  no, 
included  in  the  thermal  modelt.  cnlculutod  in  thtu  paper  h..cuu,u.  Ihe 
heat  lug  from  this  rca.tlou  was  luslgeiricaul  (Sleep.  1972). 
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Figure  Captlon» 

I 

Kuril  Islands  earthquakes.  In  this  and  subsequent  figures 
circles  are  earthquakes  and  triangles  are  seismograph 
stations.  Bathymetry  Is  In  fathons  (1  fathom  1.83 
meters). 

Velocity  model  of  Kurils  constructed  from  thermal  anomalies 
with  volcanic  arc  located  at  position  C,  Hypocenters 
and  model  are  projected  as  surface  distance  to  the  center 
of  curvature  of  the  arc. 

Station-corrected  amplitudes  on  a log2  scale  for  Kuril 
event  at  a depth  of  114  km.  Projection  lo  equl-dlstant 
azimuthal.  Lines  correspond  to  rays  at  10*  distance  and 
30*  azimuth  Increments  on  a symmetric  earth  that  have  been 
traced  through  the  plate  model  for  positions  A - D. 

Central  Aleutian  earthquakes  located  using  only  local  data 
and  plate  model  shown  In  Figure  5. 

Central  Aleutian  earthquakes  located  using  both  local  and 
telcselsmlc  data  and  plate  model  shown.  Dashed  line  Is 
hypothetical  plate  outline. 

Central  Aleutian  earthquakes  located  using  only  local  data 
and  plate  Aodel  shown. 

Composite  focal  mechanisms  for  earthquakes  at  depths 
Indicated.  Projection  Is  equal-area  azimuthal  on  lower 
hemisphere.  P and  T refer  to  prefisure  and  tension  axes, 
respectively. 

Alaska  earthquakes  and  local  seismograph  network, 

Alaska  earthquakes  located  using  plate  model  shown. 


Figure  1«10: 


Figure  l.ll! 
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Kuril  nhallow  earthquake's  demonstrating  effects  of  plate 
structure  on  focal  mechanism  determinations. 

1965  Rat  Islands  and  1964  Alaskan  malnshocks  and  plate 

A 

effects  on  focal  mechanism'  deterbinations. 
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2.  Focal  Depth  Determination  from  the  Slgtul  Character  of  Long  Period  P Waves 

Robert  B.  Herrmnn 

Introduction 

The  deTicrmination  of  focal  depths  of  earthquakes  has  lonfj  been  nn 
interesting  and  difficult  problem  in  seismology,  c;specially  for  crustctl 
events.  Some  of  the  difficulties  in  focal  depth  determination  for  fuich 
shallow  events  are  due  to  the  lack  of  sensitivity  in  the  travel  Lime 
curves  to  focal  depth  variation.  Inadequate  distribution  of  recording; 
stations,  and  a poor  signal  to  noise  ratio  for  picking  the  arrival 
times  of  the  depth  dependent  phases. 

The  use  of  short-period  seismograms  yields  the  required  time 
resolution  needed  for  precise  focal  depth  determination.  However, 
complex  scattering  effects  due  to  crustal  inhonogeneities  as  well  as 
the  earthquake  rupture  process  complicate  the  record.  Thir;  makes  the 
identification  of  the  depth  deterinit>lng  phases,  ouch  as  pP  and  nV> 
difficult.  On  the  othar  hand  long-period  seismograms  usually  have  n 
clean,  coherenr.  body  wave  signcl,  but  the  depth  phases  are  d^sl  inct 
only  for  subcrustnl  events  due  to  the  integrating  effect  of  the  long- 
period  seismograph. 

A recent  paper  by  Langston  and  llelwberger  (3  975)  sl.udic’d  the 
Ctff  cci  s of  focal  dcjitb,  focal  mechanism  and  sti-er.n  drop  upon  Uie 
recorded  long-period  telcseismic  P wave.  They  modeled  the  earth 
structure  at  the  source  and  receiver  by  a homogeneous  half space. 

Their  synthesized  teleselsmic  P v/ave  signals  contain  the  contributions 
of  just  the  direct  P w.ave,  and  the  pP  and  sP  phases.  Their  paper 
.suggests  the  possibility  of  using  the  signal  character  of  long-period 
P \7avcf.  as  a focal  dc.pfh  discvimiiiant. 


hi 

« 

The  object  of  the  present  paper  is  to  test  the  oni  l.'ibi  1 i.l  y of 
rsing  the  tcleseistaic  P-wave  signal  to  determine  focal  dcptli  using 
more  realistic  models  of  the  earth  structure  at  the  source  and 

I 

receiver  than  used  by  Langston  and  llelmberger  (19V5). 
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Theory 

The  predicted  sel&mogram  y(t)  is  determined  using  the  relation 
y(t)  “ s(t)  * R (t)  * M(t)  * R ^(t)  * I(t)  (1) 

vhere  s(t)  is  the  source  time  function  of  the  double  couple  viithout 

moment  force  system:  R (t)  Is  the  time  history  of  the  tclescismlc 

sc  ■ 

mantle  P wave  due  to  an  impulsive  source  within  the  layered  medium  iised 
for  the  crust  aiid  upper  mantle;  M(t)  Is  the  mantle  impulse  response 
for  P waves;  R (t)  Is  the  response  at  the  surface  of  the  receiver 

SC 

crust  lor  un  iticldeni.  iiopj-sivc  teleselsmlc  P wave;  1(1.)  ifJ  the 
impulse  response  of  the  seismograph  system  used;  and  * is  the 
convolution  operator. 

U.sunlly  s(t)  is  not  specif ied  when  .siudying  lelcseJsmlc.  P waves. 

Rather  n’(t),  the  time  derivative  of  s(t)  , is  used  Blnrc  the  time 

history  of  the  lar'*f.leld  body”wave  motion  is  proportional  t.o  o*(t) 

for  the  double  couple  withotit  moment  source  (Key]  i.s-Borok,  19.'30).  In 

this  study  s’(t)  is  taken  to  be  a positive  pulse,  shaped  like  an 

Isosceleo  triangle  v/Llh  a haae  of  i seconds  and  an  area  numerical] y 

equal  t:o  the  seismic  moment,  . 

R (t)  is  found  using  Hudson’s  (1969)  extension  of  the.  Haskell 
sc 

(1964)  matrix  formalism  to  the  generation  of  tclesclsmic  body  waves 
from  a point  system  of  forces  acting  within  a plane-layered  model  of 
the  upper  mantle  and  crust.  ” function  of  the  r.rufital 

model,  the  phase  velocity  of  the  telcseismlc  P wave,  the  focal  depth, 
^^10  orientation  of  the  system  of  forc’cs  at  the  source.  A point. 
douMc  rouplr  v/ltliout.  ..oment  r:o>irc'-  ur;cd  in  this  sl.udy. 
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M(t)  represents  the  impulse  response  of  the  jnnntlc  transfer 
function.  Following  Langston  and  Helmberger  (1975) > a reasonable 
approximation  for  ll(t)  in  the  epicentral  distance  range  of 
la 

M(C)  . G(4)  Q(l;)  , (?) 

where  G(A)  accounts  for  the  effect  of  geometrical  spreading  within 
the  mantle  upon  the  amplitude  of  the  P wave.  A is  the  epicentral 
distance  from  the  source  to  the  receiver.  The  values  of  G(A)  used 
in  this  paper  are  taken  from  Langston  and  Helmberger  (1975).  Q(t) 
accounts  for  the  effect  of  anelastlc  attenuation  in  the  mantle.  In 
the  frequency  domain  the  magnitude  of  this  effu  is  talccn  to  he 

|Q((o)|  *=  cxp(-toT/2Q)  (3) 

where  u is  tin.  angular  frequency,  T is  the  travel  time  of  the 
r vavii  in  Ihe  mantle,  and  Q is  the  average  Q factor  along  the  i).'n:tic.uJ  ar 
path  of  the  P wave.  Following  Langston  and  Helmberger  (1975),  we  take 
T/Q  •=  1 as  a reasonable  approximation  for  teleseismic  P waves  in  the 
epicentral  distance  range  of  A ^ 30“  to  80“ . This  epicentral  distance 
1 ange  wdj  ( fu'seii  to  avoid  signal  couiplicaLiore.  due  to  arrivals  from 
other  blanches  of  the  travel  time  curve  or  from  deep  reflections. 

In  the  present  paper  Q(t)  is  taken  to  be  non-causal.  This 
necessarily  implies  that  the  source  time  function  s(t)  , obtained  by 
fitting  tlie.  ol>served  P wave  notion  by  the  theory  of  equation  (1) , is 
not  truly  representative  of  the  history  of  the  faulting  at  the  source., 
hut  I'nther  that  «(l)  is  sucl>  tha»-  s(t)  * Q(t)  yields  the 


*.  *-J  ?1.>W-*- 
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appropriate  shape  of  the  teleseismtc  P wave,  because  of  the  natiiTfr 

of  |Q(co)|»  the  seismic  moment  determination > the  zero  frequency 

level  of  the  Fourier  transform  of  si^t),  ns  well  as  corner  froqiiency 

« • • 

information  will  not  be  affected  by  the  use  of  a non-causal  Q(t) 
filter.  The  causal  Q(t)  funct:ion  is  discussed  by  Futterman  (196Z) 
and  Stride  (1970). 

R^^(t)  gives  the  ground  notion  at  the  surface  of  the  layered 
receiver  crust-mantle  system  for  a teleseismic  mantle  . P wave 
incident  at  the  base  of  the  layered  system.  The  formalism  of  Haskell 
(1962)  is  used.  In  this  paper,  11  (l;)  is  used  to  genexate  the 

vertical  component  of  displacement  at  the  free  surface. 

The  seismograph  impulse  response  I(t)  is  .synthesized  for  a 
15--;100  iWSSli  Iv'iig  period  vortical  sj'steia  with  a peak  magni  ficati.on  of 
3000.  llnglwara's  (1958)  foriuuln  for  the  response  of  an  electroinagne*  ic. 
.‘.cir.inograph  Lu  used  together  with  the.  appropriate  system  eonntnvitf: 
given  by  Chandra  (1970). 

To  gain  an  insight  into  the  nature  of  the  focal  depth  information 

contained  in  observed  teleseismic  V wuvit  8.ii'nal,  the  contributioris 

of  R (t)  and  R (t)  to  the  observed  P wave  are  con.sidered  in 
sc  rc^ 

dettnil  for  several  earth  models. 

Tlie  earth  models  used  are  given  in  Table  1.  Model  A Is  a truncated 
version  of  the  oceanic  model  given  by  Harkrider  (1970).  Si.nco.  we  arc 
interested  in  iiiodeling  shallow  t vents,  only  the  upper  80  Uin  of 
Harkrlder’s  (1970)  oceanic  model  is  used.  Model  R Is  a clmllnvly 
truncated  version  of  the  shield  model  oh  llnrb rider  (1970).  Model  C 
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is  a siuipla  crustal  model  with  a 30  Isu  thick  crust.  Model  D has  the 


velocities  as  Model  C but  it  has  a 40  Icn  thick  crust. 


Figure  la  shows  the  fuaction  R5*'(t)  determined  by  using  the 


source  whose  parameters  arc  given  in  Table  2 under  the  event  entitled 


TEST.  Model  D of  Table  1 was  used  for  the  source  region  structure. 


The  seismograos  are  for  a ray  which  leaves  the  source  along  an  azijimth 


of  0“  with  a phase  velocity  of  16  km/sec.  Several  easily  identifiable 


phases  can  be  identified  in  the  first  ten  seconds  of  the  record.  The 


first  arrival  is  the  F wave  which  goes  downward  from  the  source  and 


becomes  the  mantle  F wave  at  the  crust-mantle  boundary.  The  second 


phase,  dencted  as  S,  leaves  the  source  as'  S and  is  converted  into 


M rorci)o.ient  of  the'  tclctielsmtc  F wave  by  Lransnilssion  at  the.  crusL- 


mntle  interface,  llie  other  phases  are  free  surface  reflections  v?hich 


are  converted  to  telcseistnic  F vmves  at  the  crust-nantle  boundary. 


it  In  uaportant  to  note  that  the  S v^avca  generated  by  the  double- 


couple source  can  make  a large  contribution  to  the  shape  of  the 


resultant  tcleseismlc  F wave. 


Figure  lb  sUo.va  the  vertical  component  of  ground  disp'J  acement  at 


the  nurliicc  of  the  receiver  crust  fox'  an  iiicident  mantle  F v;ave 


having  a phase  velocity  of  c » 16  kin/sec  and  a triangular  shaped 


pulse  with  T = 2.0  sec.  Crustal  model  D was  also  used  for  the  earth 


model.  The  phase  Identified  as  P represents  the  contribution  to  the 


vertical  displaceiuent  due  to  the  crustal  F wave  excited  by  the 


incident  mantle  F wave.  Tlie  S phase  lepresents  tlie  contribution 


due  to  the  mantle  F wave  to  crustal  S vrave  conversion  at  the  cru.'it- 
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mantle  interface.  Ihe  PPP  arival  represents  tlie  contr.i.bution  ilue  to 
the  first  reverberation  within  the  crust  !.ur.h  that  the  wave,  travels 
as  P on  call  three  segments  of  the  ray  path.  The  next  arrival  cm 

I » 

the  seisDOgram  trace  is  due  to  the  siiaull  aneous  arr.lval  of  the  throe 
rays  SPP,  PSP,  and  PPS.  It  is  seen  that  these  multiple  ref locti.fins 
arc  large  enough  to  affect  the  character  of  the  resultant  surface 
displacement  up  to  twenty  seconds  after  the  direct  P arrival. 

The  relative  time  delays  between  the  arrivals. is  a function  of  the 
velocity  model  ^rithio  tue  crust.  The  partition  of  amplitudes  l)otween 
the  P and  S contributions  is  basically  a inncLlon  of  l)ie  veJocity 
and  density  contrasts  at  the  crust-nantle  boundary.  'Jlie  crustal  model 
o hore  has  a very  piononnced  miterial  c.oaLv-'if.t  nl  the  f / unt-mant  li- 
hemuinry.  Ihiir.  t.ho  1’  to  S ami  S to  P convcrcions  at.  the  crusl.- 
KK’.ntlo  boimdnr)’  r'ay  he  more  efficient  for  this  s:hn))lc  earth  model  than 
lor  ;i  i iu'f  di'J  n Ll«.d  I’.'irth  model. 

Figure  Ic  gives  Mie  composite,  effect  of  both  the  noiircc.  and 
receiver  crusts,  (obtained  by  convolving  the  rricord  sections  of  J'ignrcn 
la  ...id  lb.  The  t.irsL  ten  to  tv.el',  e .scc.oiuIk  of  fiignaJ  1 ol  1 cjwi  ng  the 
P .irrlvat  of  i.ignie  ic  are  cpiil.,'  siuflai  to  tho..f.  of  f'ig,iie  la 
_ha  receiver  crust  does  not  pertuib  the  tele.seisniLc  P \.ave  signal 
TTUch  in  this  Line  frame.  However,  the  receiver  crust  does  affect  the 
n.Vc.uTc  of  the  incident  telcseir.nlc  P wave  signal  .at  later  times. 

Pignre  2 stiows  the  effect  (»f  .focal  depth  find  dlllcri'iil.  c.nislal 

modeling  upon  th.'  oh  •.etved  vertlcaJ  displacement  oi  the  P wave; 

siin.n.  '.riie  traces  ;hovai  .n  e tin  functi.o.i  s(t)  * It  (I;)  <<■  U (t) 

* * fi  C I c 


hovai  .n  e tin  functi.o.i 
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for  the  same  source  and  receiver  conficuration  used  in  Figure  1. 

Figure  2a  uses  Kodel  C for  both  the  source  and  receiver  crusts. 

Figure  2b  used  Model  D,  and  Figurfe*  2c  lias  Model  D for  the  source  and 
receiver  crusts.  The  synthetic,  seismograms  were  gennirated  for  focal 
depths  of  10  and  20  kilometres.  The  effect  of  focal  depth  upon  the 
signal  character  is  quite  profound.  In  Figure  1,  it  was  seen  that 
the  prlae  contributoza  to  the  teXeseisaic  F wave  were  the  P,  pP 
and  aP  phases.  The  time  delays  of  the  surface  reflected  pP  and 
sP  phases  relative  to  the  direct  P phase  are  functions  of  the  focal 
depth.  For  shallow  events,  the  time  delays  are  small.  For  the  • 
particular  choice  of  the  source  orientation  used  here,  the  pP  nnd 
sP  idiases  luve  tlie  sane  polarity  so  that  they  interfere  constructivcJy 
for  the  10  km  evfnt  to  produce  a broad  arrival  following  the  direct 
P contribution.  As  the  focal  depth  increases,  the  phase  separation 
lucreuses  mid  the  pP  and  sP  phases  appear  as  disLiuct  arrivals. 

The  effocc  of  different  continental  mist  models  il;:  not  as 
profound  ns  the  .focal  depth  ef  fect.  Comparing  a s.lriplG  cru.st  with  a 
crusta.l  thickness  of  30  lire.  Figure  2a,  with  one  with  a 40  km  thick 
crust.  Figure  2b,  no  major  difference  is  seen  other  than  a slight 
change  in  the  character  of  the  latter  part  of  the  trace.  Contrasting 
Figure  2c  \rfth  Figure  2b,  the  effect  of  a more  detailed  crustal  model 
is  a smoothing  of  waveform  shape  for  arrivals  after  the  direct  P 
contribution  doe  to  the  presence  of  many  more  crustal  multiples  in 


the.  record  section. 
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Observations 

In  Equation  1,  we  note  that  many  variables  affect  the  shape  of  the 
ohserved  seismogram.  These  factors  are  the  orientation  of  the  earth- 
quake focal  mechanism^  the  source  time  function^  the  crustal  structures 
at  the  source  and  the  receiver,  and  the  focal  depth.  TInnally,  a rough 
estimate  of  the  crustal  structure  can  suffice  in  modeling  the  P arrival 
since  we  have  seer  from  Figure  2 that  the  crustal  structure  does  not 
affect  the  gross  detail  of  the  predicted  ground  motion  very  much.  On 
the  other  hand,  the  orientation  of  the  focal  mechanism,  the  shape  of 
(he  source  time  function,  and  the  fornl  depth  are  usually  not  very  well 
known.  Various  search  techniques  or  inversion  techniques  have  been 
used  to  dctcrmice  these  parameters  by  finding  the  bast  fit  between  the 
predicted  and  observed  seismograms  at  a number  of  stations.  Some: 
useful  techniques  have  been  outlined  by  Mcllman  al.  (197.'>)  and 
bangston  (19V.5). 

Since  the  object  of  this  paper  is  to  demonstrate  the  potential  oi 
using  long  period  P v?aves  for  focal  depth  resolution,  we  reriuiro  that 
the  focal  inrchanlr.ra  ami  source  time  function  as  well  as  focal  depth  he 
dcJcriiilned  previously.  This  detcrnination  \/as  made  by  finding  the 
combination  of  these  parameters  which  gave  the  best  fit  to  the  P wave 
signal  at  fifteen  selected  \WSSN  stations.  The  parameters  of  the  tx^o 
earthquakes  so  considered  are  giveii  m Table  2.  Tire  source  parfimeters 
given  are  the  orientations  of  the  pressure  and  tension  axes,  the  base 
length  of  the  isosceles  triangle  shaped  source  time  function,  the 
seismic  moment  and  the  focal  depth.  All  the  factors  given  in  Equation 
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(1)  are  Included  in  the  synthesis  ot  the  seismograms. 

The  first  event  considered  is  the  Alma-Ata  event  of  June  !>,  19/0. 
Tills  event  occurred  at  04.53:06  with  can  epicenter  42.5“  H and  78.8''  E. 
Figure  3 presents  the  observed  vertical  component  long-period  I'-uav.^ 
signal  at  CHG,  a distance  of  29.2“  from  the  so-ircc  along  an  azimuth 
of  13S“.  Also  presented  in  this  figure  are  the  synthesized  long-period 
P wave  signals  for  various  focal  depths  using  earth  model  B for  both 
the  source  and  receiver  crusts.  At  this  eplcentral  distance,  the  P " 
wave  phase  velocity  for  focal  depths  between  0-50  l;m  is  insensitive 
to  focal  depth  and  a value  of  11.02  km/sec  was  used  (Herviii  et^  a^l,  :l9f.8). 
The  focal  depth  range  of  the  synthetic  Beisnogtams  is  between  1 and 
5.1  km.  All  other  .source  para’seters  are  held  const nnt.  'flic  .scaling 
bar  represents  1 ca  of  dispiacer:eat  on  h 1.5-100  VIW.58N  seismograph  with 
a peak  nagnific ation  of  3000. 

‘I'l’.e  ofieef  of  focal  depth  is.  exhibited  in  sicveraJ.  v^.tys.  A;-,  the 
for.nl  depth  .incre.nses  within  tlie  crust,  some  character  is  introduced 
into  the  first  part  of  the  s-iigual  due  to  the  incre.nsLng  pP-P  and  .sP-P 
!.uue  <!elsrys.  As  (Ue.  focal  depti-,  .iucieas'.es  so  that  the  source  lies,  in 
ih<.  npiisv  III  iiClr,  il"u  diri'ct.  i wav<;  umplitudc  d'  Ci'i’.i  ; ber.iiisr-  the 
direct  P wave  contribution  leaves  the  source  ne.nr  a nodal  plane  (in 
other  words,  since  the  phase  velocity  does  not  change,  the  takeoff 
angles  are  functions  of  the  .source  layer  velocity  and  hence  focal 
d pth).  A visual  comparison  of  the  olseirved  (JIG  record  that  Lhe. 
synthetic  records  indicates  that  the  focal  depth  of  the  event  v;onld 
i>e  I'l.iekcliul  Ih  UU'I'U  H 1 l.n  ;md  II  21  km. 
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To  obtain  a wore  quantitative  choice  of  the  beat  fit  tiny  focal 
depth,  the  following  procedure  xms  used.  A pair  of  the  twelve  traces 
Cwen  in  Figure  3 were  cross-correlated  over  a time  .shift  of  .1/, 
s^econds  from  the  beginnitTg  of  the  trac;.e.  hTien  the  time  shift  giv.mg 
the  maximum  cross-correlation  was  found,  a correlation  coefficient 
between  the  shifted  traces  was  determined;  a correlation  coefficient 
equal  to  1.0  indicates  excellent  fit  between  the  two  traces.  The  cross- 
correlation  was  used  because  iu  actuality  there  is  always  .some  doubt 
as  to  the  exact  arrival  time  of  the'  P pulse  and  in  tT.e  inifi.a.l on  of 
the  digitizing  process  of  the  observed  signal. 

Table  3 presents  the  correlation  coefficients  obtained  by  correlatiug 
the  first:  37  r;ccouds  of  the  trcces  sh.ncu  Since  Cite  rwaj.rix  of  corrcl- 
nKon  coefficMuits  is  .symmetric  about  the  diagonal,  only  half  of  the. 

VK.trlx  Ir  given  in  Table  3.  ConsIderJog  the  correlations  hetween  the 
^ -Imu;  synthetic  ( races,  ncc  rhat  Llievc  is  a good  ecu  relat  :i.on 
{’.reatc^r  than  0.900,  for  focal  depths  crLthin  ilO  km  of  the  chosen 
focal  depth.  On  the  other  hand,  the  matrix  indicates  a distinct  non- 
corret.  Lion  he  I ween  any  two  traces  whose  focal  depths  are  separated 
h-  least  15  Ku.  Compariug  ti.  correlations  between  the  observed 
Clt.  trace  and  those  for  the  various  focal  depths,  it  in  seen  that  the 
best  fit  occurs  in  the  neighborhood  of  H = 16  Ion,'  with  anything 
deeper  being  positively  rejected.  In  order  to  set  confidence  limits 
on  the  focal  uepth.,  a gcnctralized  inversion  technique  wc;uld  have  to 
be  used  CMellman  et  19/5).  Ibis  event  has  been  assigned  focal 

I ’ 

cb-I.rir;  of  il  20  kn  by  MlA/v  and  K =■  7M  km  by  the  ISC. 
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Thu  problem  of  focal  depth  determination  for  eventr.  locaLed 
within  an  oceanic  crust  was  also  considered.  'Hie  Aleutian  Island  i 

event  of  August  7,  1966  was  used.  The  focal  neehanisro  paramete.ra  listed 

I 

in  Table  1 are  taken  from  Stauder  (1968).  Tlie  origin  time  Is  07:11;04 
and  the  epicenter  is  located  at  50.6“  N and  171.3"  U.  A truncated 
version  of  the  oceanic  model  of  llarkrider  (1970),  Model  A of  Table 
was  used  for  the  source  crust.  The  continental  crust  model.  Model  B 
was  used  for  the  receiver  crust  at  KEV.  The  seismic  moment  and  source 
pulse  shape  were  determined  in  this  study.  KEV  is  59.2"  from  the 
source  along  rui  iizlmuth  of  352.7".  A phase  velocity  of  16.092  km/f;ec 
was  used.  The  observed  KEV  P wave  signal  is  given  in  Figure  4 
together  vrjth  synthetic  records  for  various  focal  depths.  The  scaling 
bar  represents  2.5  cm  displace.iirtnt  on  :i  15-100  .1  ong -imri t)cl  f!e.ir.mogram 
with  peak  magnification  of  3000. 

A conmarir.on  of:  Figure  4 x/ith  Figure  3 point.*;  up  n proioimd 
difference  in  the  character  of  the  theoretical  long  period  P wave 
signal.  Tuis  difference  is  the  almost  consLant  period  reverberation 
following  the  iuii  ial  F wave  .u.el  lon  for  the  i venl  occ  ui’T.Ing,  ;Iu  the. 

| , oceuaic  ii’odel.  This  reverheratioa  is  due  to  the  prcsem-ci  of  the  \’.ite.r  ] 

layer,  v;i.th  the.  period  of  the  reverberation  being  a function  of  I he. 
thickness  of  the  water  layer. 

The  first  30  seconds  of  the  seismograms  in  Figure  4 were  correlated 
using  the.  !:.'uie  procedure  as  for  the  CHO  record.  'Jhe.  fo.llov/iiig  u-jura 
correlation  coefficients  were  determined  between  KF.V  and  the  various 
focal  de.ptlis:  0.526,  0.897,  0.961,  0.967,  0.974,  0.903,  0.821,  and 


0.A84  for  focal  depths  of  6.5,  10,  I5,  .19.5,  24.5,  27.5,  30,  and  35  lun, 
respectively.  Even  though  the  H “ 6.5  km  correlation  coefficient  .Is 
quite  good,  the  theoretical  P wave  break  is  not  as  sharp  as  that 
observed  on  the  KEV  record.  On  the  basis  of  these  correlation 
coefficients  and  a visual  comparison  of  the  traces,  a focal  deyjth  in 
the  range  of  15  - 24.5  km  would  explain  the  KEV  seismogram.  Performini 
this  analysis  x>rith  a number  of  observed  seismograms  would  provide  n 
tighter  control  on  the  focal  depth.  For  comparison,  the  focal  depth 
of  this  event  was  given  as  H “ 39  km  by  the  USCGS  and  H “ 29  km  by 

the  ISC. 
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Conclusions 

flic  po£isibility  of  using  the  signal  character  of  long-period  P 
\7aves  has  been  demonstrated,  but  something  must  be  said  about  the; 
applicability  cf  this  technique.  The  present  study  indicates  that  it 
is  not  unreasonable  to  obtain  a resolution  of  about  fifteen  kilometers 
In  focal  depth.  This  means  that  an  event  could  be  classified  as  having 
occurred  in  the  upper  or  In  the  lower  crust. 

There  are  some  restrictions  in  obtaining  this  resolution.  First 
and  foremost  is  tba  quality  of  the  data  available.  The  events  considered 
here  have  a relatively  high  signal  to  noise  ratio.  The  six-second 
microseism  level  may  be  on  the  order  of  several  tenths  of  a centimeter  •• 
on  a typical  15-100  V'if'SSM  ceisnogr.'im  \d.th  a magiilficaUion  of  3000.. 
llct\.’ecm  the  microseioca  noise  leve3_  and  the  width  of  the  seismogram 
trace,  it  v/ould  l»e  cUtficult  to  apply  this  technique  to  events  v/Ith 
i.c-Lsmic  i.'orjpnts  loss  than  one- tenth  th.'iU  of  the  dune  5,  V)/{)  event. 

v.liicU  luid  an  iL  ■ 6.3  and  an  m ” 6.0. 

I*  ■ i> 

The  P - \;ave  detectability  nay  be  improved  by  using  a high-g.ain 
long-period  seismograph  system  designed  to  reject  the  mictoselsm  noise, 
hut  tlie  lesult  Juay  he  a greater  lack  of  resoluLion  in  the  depth 
resolution  due  to  the  fact  that  a longer  period  seismograph  is  being 
used.  This  aspect  of  the  problem  \-rill  be  covered  in  a future  study. 

Since  the  emphasis  of  this  study  was  toward  focal  depth  resolution 
for  small  events,  the  problem  of  a complex  three  dimensional  source 
V7as  not  considered.  However,  the  v/ork  of  Heilman  ct  al  (197.5)  and 
T.ane.ston  (197.5)  iud i rate.'?  that  inf  orir.ntiuu  coiiceiming  depth  .and  rupt  uring 
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can  be  extracted  irom  the  long-period  telcseismic  V wave. 

The  present  study  may  also  be  described  as  serai-quantitative,  it 
Xfould  be  well  to  study  this  problem  from  the  viewpoint  of  generalized 
inverse  theory  in  order  to  be  able  to  set  some  confidence  limits  on 
the  focal  depths  determined. 
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TABLE  3 

CbrrelatloQ  matrix  for  CHG  record  05  JUN  VO 


CHG  Z 

1 km 

6 kra 

11  km 

16  km‘ 

21  kra 

31  km 

41  km 

51  km 

CHG  Z 

1.000 

.868 

.918 

.934 

.958 

.894 

. 286 

.230 

.248 

1 kra 

1.000 

.981 

.976 

.916 

-.087 

-.035 

.022 

.146 

6 km 

1.000 

.995 

.983 

.358 

-.186 

.040 

.146 

11  kn 

■ • 

1.000 

.989 

.738 

-.188 

.089 

.135 

16  km 

1.000 

.961 

.102 

.184 

.212 

21  km 

1.000 

.938 

.152 

.300 

31  km 

• 

1.000 

.944 

.342 

A1  Ian 

1. 000 

.997 

51  km  . 

1.000 
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Figure  Captions 


Figure  2.1. (a)  Mantle  P-wave  due  to  source  crust  interaction  1 
(b)  Vertical  displacement  at  free  surface  of  receiver 
crust  due  to  an  incident  nantle  P-wave  1 (c).  Convolution 
of  source  crust  and  receiver  crust  responses. 


Figure  2.2.  p!f feet  of  focal' depth- and  crustal  structure 
upon  the  resultant  vertical  component  ground 
displacements.  (a)  earth  model  C i (b)  earth  model  D 1 


(c)  earth  model  B. 


Figure  2.3.  Observed  P-wave  signal  at  CHG  for  the  June  5,  I970 
event  compared  to  synthetics  for  various  focal  depths. 


The  scaling  bar  v represents  1 cm  displacement  on  a 


3K  15-100  V7WSSN  seismogram. 


Figure  2.4.  Observed  P-wave  signai.  at  KEV  for  the  August  7,  I966 
event  compared  to  sjnthetics  for  various  focal  depths. 

The  scaling  bar  represents  2.5  cm  .diuplacoment  on  a 


3K  1 5-1 00  vr.i’SSN  seismogram. 


I 


1?  1 7.1 
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3 . Source  Characteristics  fr'~*r>  I.ong  Period  P Wave  Observations 

R.  B.  Herrmann,  and  T.  Hewitt 

In  the  previous  semi-anriual  report,  we  presented  a method  determining 
earthquake  source  parameters  from  a study  of  the  long->period  teleseismic 
P-wavc  character.  The  theoretical  dlsucssion  in  that  report  has  been 
updated  by  using  the  corrected  reflection  and  transmission  coefficients 
for  seismic  rays  in  the  source  region  (Bache  et  al..  1975b).,  As  a further 
extension  of  the  theory  given  in  the  previous  semi-annual  technical  report, 
we  used  Hudson's  (1969)  extension  of  the  Haskell  (1962,  1964)  formalism 
to  produce  synthetic  seismograms  for  teleseismic  P waves  for  source 
crust  and  receiver  crust  models  which  are  more  complicated  than  the 
single  la3rered  constant  velocity  crust  overlying  a constant  velocity 
mantle,  used  previously. 

The  technique  used  in  determining  the  source  parameters  of  an 
earthquake  is  quite  syjnple.  Long-period  teleseismic  P waves  fromstations 
located  uniformly  about  the  epicenter  are  compared  with  theoretical  P 
waver,  generated  for  a particular  combination  of  strike,  dip,  slip,  focal 
depth  and  source  time  function.  The  parameter  space  of  source  parameters 
is  searched  until  a combination  is  found  which  provides  the  best  fit 
to  both  the  amplitude  and  the  shape  of  the  observed  P-wave  signals. 

Twelve  earthquakes  have  been  studied.  Table  1 lists  the  events, 
e.g.,  giving  the  date,  origin  time,  epicentral  coordinates,  focal  depth 
and  magnitudes.  The  focal  depths  and  magnitudes  listed  were  determined 
by  the  CGS  or  NOAA,  the  ISC  and  this  study.  No  depths  were  determined 
for  events  1 and  2,  Aleutian  trench,  because  the  search  technique  is  set 
up  for  continental  events.  Events  4,  5,  and  12  are  too  deep  to  be  processed. 
The  focal  depths  determined  from  the  teleseismic  P-wave  analysis  differ 
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somewhat  from  those  given  by  the  CGS  or  NOAA  and  the  ISC.  Our  depth 

estimate  of  the  22  May  1971  Turkish  event  is  the  only  one  which  differs 

substantially  from  the  estimates  of  the  other  two  sources.  We  also 

• • 

determined  the  magnitude  of  the  events  together  with  standard  deviations. 
The  number  in  the  parentheses  indicates  the  number  of  data  used  in  our 
determinations.  In  general,  our  magnitude  determinations  are  somewhat 
less  than  those  from  the  other  sources,  but  the  other  determinations 
usually  fall  within  our  standard  deviations. 

Tabic  II  presents  the  focal  mechanism  parameters  of  the  events 
JisTcd  in  Tabic  T.  Of  the  12  events  listed,  we  determined  the.  focal 
mechanisms  of  seven  events.  The  P-wave  first  motion  data  and  S-wave 
polarization  data  are  presented  in  Figure  1,  together  with  the  focal 
mechanisms  determined  from  the  long  period  ?-wave  signal  analysis.  In 
these  figures,  the  O's  represent  compressional  arrivals  while  the  deltas 
represent  dilatational  arrivals.  The  X*s  represent  arrivals  of  near 
nodal-plane  quality.  The  P and  T symbols  represent  the  pressure  and 
tension  axes.  The  event  of  12  Nov  1972  was  a deep  event  and  readings 
were  available  from  the  pP  arrival;  the  minus's  represent  dilatations 
obtained  from  the  pP  arrivals  after  correction  for  the  phase  reversal 
due  to  the  free  surface  reflection. 

The  focal  mechanism  of  the  lA  Oct  1968  Meeker Ing,  Australia  agrees 
relatively  well  with  the  observed  surface  faulting.  The  surface  faulting 
extended  a distance  of  30  Ion  in  an  arcuate  trend,  striking  southeast 
at  the  southern  end  of  the  fault  and  northeast  at  the  northern  end. 

The  overall  trend  is  north-south.  The  motion  on  th»*  fault  was  a thrust 
along  a plane  dipping  30®  - A9*  at  the  surface  (Gordon,  1971). 

Banghar  (1974)  published  focal  mechanism  solutions  for  the  11  Feb 
1969,  18  Jul  1969  and  the  05  Jun  1970,  events.  Ills  solution  for  the 
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11  Feb  1969  event  has  one  nodal  plane  dipping  70“  to  the  SW 
with  a strike  of  38“  and  the  other  nodal  plane  dipping  58*  to  the  SW 
with  a strike  of  321“.  There  is  no  similarity  whatsoever  between  our 
solution  and  his.  The  solutions  for  the  05  Jun  1970  events  are  quite 
similar.  One  of  his  nodal  planes  dips  50*  to  the  NW  with  a strike  of 
68*  while  the  other  nodal  plane  dips  50*  to  the  SE  with  a strike  of  26*. 

This  is  very  similar  to  our  solution,  the  difference  being  a 20*  difference 
in  the  strike  of  the  nodal  planes.  This  discrepancy  is  not  major  since 
we  found  the  P-wave  signal  to  be  relatively  insensitive  to  a rotation 
of  the  focal  mechanism  strikes  since,  amplitnden  of  the  ? waves  and  SV 
waves  are  fairly  insensitive  to  a rotation  of  the  solution  for  the  rays 
making  up  the  telesclsmlc  P wave. 

lianghar's  (1974)  solution  for  the  18  Jul  1969  northeastern  China 
event  has  one  nodal  plane  striking  along  an  azimuth  of  109*  and  dipping 
70*  to  the  south,  and  the  second  nodal  plane  striking  along  an  azimuth 
of  22“  and  dipping  80*  to  the  west.  This  solution  is  quite  rirailar  to 
ours.  The  P-wave  signal  comparison  technique  began  to  break  down  for 
this  event  since  this  event  was  large  enough  that  the  assumption  of  a 
point  d-*3location  source  began  to  break  down. 

Besides  determining  the  focal  depth  and  nodal  plane  ox'ientatlons, 
we  were  able  to  estimate  the  shape  of  the  body  wave  source  time  function. 

Since  we  were  using  teleselsmic  P waves  from  a 15  — 100  WWSSN  seismograph 
and  because  of  the  effect  of  the  earth’s  Q filter,  we  lack  sensitivity 
in  the  fine  structure  of  the  source  time  function  which  can  be  obtained 
only  from  the  high  frequency  content  of  the  signal.  To  try  to  resolve 

i 

the  nature  of  the  high  frequency  content  of  the  signal,  we  obtained 
the  shapes  of  the  rectangular  and  isosceles  triangular  source  pulses 
which  gave.  the.  host  fit  to  the  d.ita. 
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Table  III  presents  the  results  of  these  detenuinat jons.  In  this 
table,  TMAX  is  the  number  of  seconds  of  P-wave  signal  used  In  the  deter- 
mination of  the  best  fit;  N Is  the  number  of  stations  making  up  the 
P-wave  data  set;  H Is  the  focal  depth  ^an  asterisk  Indicates  that  the  focal 
depth  was  not  determined  as  part  of  the  search  technique);  DTI,  DT2  and 
DT3  represent  the  time  lengths  of  the  positive,  zero  and  negative  slopes 
of  a positive  trapezoidal  source  time  function.  The  case  DTI  ■>  DT3  and 
DT2  0 represents  an  Isosceles  triangular-shaped  source  time  function 
while  the  case  DTI  - DT3  - 0 and  DT2  ^ 0 represents  a rectangular 
shaped  source  pulse.  is  the  seismic  moment  In  dyi\e-cm  determined 

for  each  source  time  function,  r Is  the  correlation  coefficient  between 
the  observed  and  predicted  P-wave  signals  for  the  first  TMAX  seconds 
for  all  N data  sets. 

In  general,  the  seismic  moment  determlnstions  are  Independent  of 
the  source  time  function  assumed.  This  is  true  so  long  as  the  corner 
frequency  associated  with  the  source  time  function  is  greater  than  the 
minimum  frequency  resolvable  in  the  P-wave  signal.  Since  the  focal  depth 
determination  depends  on  the  Interference  of  the  primary  and  reflected 
body  wave  contributions  to  the  tclesclsmlc  P wave,  the  focal  depth 
determinations  is  dependent  upon  the  sliapo  of  the  input  signal  as  well, 
as  the  time  delays  due  to  focal  depth.  Hence,  the  focal  depth  estimate 
will  be  somewhat  dependent  upon  the  shape  of  the  source  time  pulse. 

Figure  2 presents  a plot  of  the  seismic  moments,  , versus 

the  surface  wave  magnitudes,  , detemdued  in  this  study.  The  niimbers 

refer  to  the  event  number  in  Tables  I and  III.  The  error  bars  give  the 

standard  deviations  of  the  M determinations.  While  standard  deviations 

8 

of  the  M determinations  were  not  obtained,  they  are  probably  on  the  order 
o 

of  a factor  of  2.  The  dashed  line  Is  an  extrapolation  of  seismic  moment 
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versus  M data  for  the  central  United  Stated  (Street  et  al. » 1975) . 

8 

The  light  solid  line  Is  the  versus  relation  based  on  Akl^s 

(1972)  model  of  source  spectrum  scaling.  Akl’s  model  provides 

a good  fit  to  most  of  the  data.  Exceptions  are  the  three  Aleutian 

Trench  events  of  01  Oct  1965,  02  Jun  1966  and  07  Aug  1966,  which  have 

consistently  smaller  values  than  indicated  by  the  versus 

M trend.  The  next  section  will  show  why  we  attribute  this  discrepancy 
s 

to  M rather  than  M . 
s o 

As  an  attempt  at  discerning  the  high  frequency  character  of  the 
source  spectrum  of  these  events,  we  considered  the  relationship  between 
the  spectral  amplitude  level  at  1 Hz  and  m^  for  each  of  the  two  source 
time  functions  from  each  event.  The  amplitude  spectrum  of  a rectangular 
pulse  of  duration  DT2  seconds  has  a corner  frequency  at  « ir/DT2  Hz, 

which  separates  the  f°  and  f"^  trends.  For  an  isosceles  triangular 
shaped  source  time  function  with  a base  of  2*  DTI  seconds,  the  spectrum 
lias  . corner  frequency  at  f(,2  " » ''^ich  separates  the  f®  and 

trends  of  the  spectrum.  Obtaining  the  comer  frequency  and  seismic 
moment  from  Table  111  for  each  source  time  function,  and  knowing  the  high 
frequency  spectral  trend,  the  spectral  amplitude  at  a frequency  of  1 Hz 
is  easily  determined.  Table  IV  sumnarlzes  these  computations.  In  this 

table,  Tqj^  = 1/^01  ^02  “ ^^^02  . 

Figure  3 shows  the  spectral  level  at  1 Hz  as  a function  of  m^  . 

In  this  figure  the  numbers  refer  to  the  events  listed  in  Table  IV.  The 
open  triangles  indicate  the  spectral  level  appropriate  to  a rectangular 
source  palse  while  the  solid  circles  indicate  the  spectral  levels  obtained 
rom  the  triangular  shaped  source  pulse.  The  error  limits  are  the  standard 
deviations  on  tho  values  listed  in  Table  I.  The  spectral  levels 

are  probably  luiown  to  within  j factor  of  2,  due  to  a similar  uncertainty 


JA 


in  the  seismic  moment  values. 
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The  dashed  line  represents  the  extrapolated  trend  of  the  1 Hz 

spectral  level  as  a function  of  aa  determined  by  Street  al.  (1975). 

It  is  interesting  to  note  that  the  best  fit  to  this  trend  is  from  an 

assumed  triangular  shaped  source  pulse  rather  than  a rectangular  shaped 

pulse.  From  this  simple  comparison,  vte  conclude  that  the  source  spectra 

_2 

near  1 Hz  approximately  have  an  f trend  for  most  of  these  events. 

A possible  exception  might  be  events  7 and  11  which  may  have  an  f”^’^ 

spectral  trend  near  1 Hz.  Events  4 and  5 are  deep  focus  events  which 

_2 

indicate  that  a falloff  faster  tlum  f might  be  more  appropriate, 
ihin  indicates  either  that  the  corner  frequency  of  these  two  events  are 
probably  not  very  wall  kno%m  or  that  a source  time  function  is  somewhat 
smoother  than  a triangular  shaped  pulse.  Events  1,  2,  3,  the  Aleutian 
Trench  events,  fit  the  trend  nicely  i>ith  a triangular  shaped  pulse. 

This  indicates  that  the  value  is  fairly  good  and  that  tht;  discrepancy 

in  Figure  2 Is  due.  to  N . 

8 

In  couclurlon,  the  following  observations  are  made  concerning  the 
use  of  long-period  teleselsmlc  P waves  for  source  characterizations. 

a)  The  capability  of  focal  mechanism  determination  is  quite  good 
since  the  observed  long  period  P-waves  train  is  composed  of  samples 

of  both  I’  and  SV  waves  generated  by  the  source. 

b)  Focal  depth  discrimination  is  good.  This  topic  is  convered 
elsewhere  in  this  report. 

c)  The  shape  of  the  source  spectrum,  estimated  from  the  solution 
of  the  P-wave  signal  problem  agrees  well  with  spectral  scaling  relations. 
From  the  source  spectrum  scaling  of  1 Hz  with  m^  , we  concluded  that 
the  triangular  shaped  source  pulse  Is  a better  solution  than  n rectangular 
puliio.  Note  that  we  had  to  use  lu^  to  be  able  to  make,  this  dlsl  liicl  Imi. 
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Recent  programming  changes  have  made  the  method  of  determining  source 
characteristics  simpler  to  use,  and  more  economical  of  computer  time. 

The  optimization  of  computer  codes  has  been  achieved  by  using  the  fact 
that  (all  other  source  parameters  fixed)  the  seismograms,  and  hence  their 
correlation  coefficients  and  seismic  moments,  are  linear  functions  of  the 
amplitudes  of  the  four  source  rays  (2P  and  2SV) , which  are  determined 
by  focal  mechanism.  Because  of  these  changes  computing  efficiency  has 
been  Improved  by  a factor  of  five  to  100  times  for  various  parameters. 

This  Improvement  in  computing  efficiency  has  allowed  most  of  the  search 
for  source  parameters  to  be  done  automatically,  and  hence  total  time 
to  process  an  event  is  now  reduced  to  only  two  days.  Instead  of  two 
weeks. 

I'uture  atudle.s  should  encompass  the  following; 

a)  Estimation  of  the  lower  magnitude  threshold  for  the  application 
of  this  technique; 

b)  Study  ol  the  1 Hz  spectral  character  from  short,  period  P waves 
to  test  the  validity  of  our  conclusions  of  the  shape  of  the  source  pulse 
for  these  events; 

c)  Investigation  of  the  cause  of  the  relatively  low  value  for 

the  Aleutian  Trench  events. 
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TABLE  II 

FOCAL  MECHANISM  PARAMETERS 


Event 

Date 

I. 

Axis 

P Axis 

Ref. 

Trend 

Plunge 

Trend 

Plunge 

1. 

01 

Oct 

65 

182 

12 

66 

65 

Stauder  . 

(1?68) 

2. 

02 

Jun 

66 

223 

5 

353 

82 

Staudcr 

(1968) 

3. 

02 

Aug 

66 

331 

64 

200 

18 

Stauder 

(1968) 

A. 

15 

Feb 

67 

257 

10 

77 

80 

5. 

09 

Sep 

67 

263 

20 

64 

69 

6. 

lA 

Oc.f 

6S 

27 

59 

271 

15 

7. 

11 

Feb 

69 

. 235 

69 

342 

7 

Stauder 

(1975) 

8. 

18 

Jul 

69 

186 

12 

91 

22 

Stander 

(1975) 

9. 

05 

•Uin 

70 

2A2 

71 

338 

2 

10. 

22 

May 

71 

100 

4 

190 

4 

11. 

03 

Sep 

72 

175 

83 

39 

5 

12. 

12 

Nov 

72 

236 

10 

346 

62 

c 

I 

■i 

. :J 

I 


TABLE  111 
GOODNESS  OF  FIT 


Date 


Ox  Oct  03 

02  Jun  66 
07  Aug  66 
15  Feb  67 

09  Sep  67 

14  Oct  68 

11  Fob  69 
18  Jul  69 
05  Jun  70 
22  May  71 

03  Sep  72 

12  Nov  72 


DTI  DT2 


o 

00 

E 

:ij 

6.0 

E 

24 

1.5 

E 

26 

4.2 

E 

26 

4.4 

E 

26 

3.f> 

E 

26 

1.9 

E 

26 

6.0 

E 

25 

4.5 

E 

25 

.2  E 25 


2.2 

E 

25 

2.0 

E 

26 

2.0 

E 

26 

6.0  E 

26 

6.0 

E 

26 

1.2 

E 

26 

1.2 

E 

26 

1.7 

E 

25 

1.7 

E 

25 

5.0 

E 

25 

5.0 

E 

25 
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TABLE  IV 

S (1  Hz)  SPECTRAL  ESTIMATES 

1 « 

Triangular  Impulse  Rectarpqilar  Impulse 


Event 

Date 

M 

o 

dyne- 

•cm 

T 

02 

sec. 

S (1  Hz) 

M 

o 

dyne- 

-cm 

T 

02 

sec. 

S (1  Hz) 

1. 

01 

Oct 

65 

8.0 

E 

25 

7.5 

1.4 

E 

24 

2. 

02 

Jun 

66 

6.0 

E 

24 

4.1 

3.6 

E 

23 

3. 

07 

Aug 

66 

1.5 

E 

26 

9.4 

1.7 

E 

24 

4. 

15 

Feb 

67 

4.2 

E 

26 

11.3 

3.3 

E 

24 

4.4  E 

26 

17.3 

2.5 

E 25 

5. 

S' 

Sep 

67 

i.'d 

K 

26 

6.3 

4.5 

h 

>4 

1.9  )i 

26 

11.0 

1.7 

K 25 

6. 

14 

Oct 

68 

5.0 

E 

25 

9.4 

5.7 

E 

23 

4.5  E 

25 

11.8 

3.8 

E 24 

7. 

11 

Feb 

69 

2.2 

i:: 

25 

9.7 

2.3 

E 

23 

2.2  E 

25 

14.9 

1.5 

E 24 

8. 

18 

69 

2.0 

E 

2,6 

12.9 

!.? 

E 

26 

2.0  F. 

26 

18.8 

1.3 

E 25 

9. 

05 

Jun 

70 

6.0 

E 

25 

11.6 

4.4 

E 

23 

6.0  E 

26 

18.8 

3.2 

E 24 

10. 

22 

May 

71 

1.2 

E 

26 

9.4 

1.4 

E 

24 

1.2  E 

26 

12.6 

9.5 

E 24 

11. 

03 

Sep 

72 

1.7 

E 

25 

6.3 

4.3 

E 

23 

1.7  E 

25 

9.4 

1.8 

E 24 

12. 

12 

Nov 

72 

6.0 

E 

25 

5.7 

1.9 

24 

6.0  E 

25 

9.4 

6.4 

E 24 
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Figure  Caplions 

Figure  3.1:  Focal  mechanism  solutions  made  as  part  of  this  study.  The 

source  parameters  are  given  in  Tables  I and  II. 

Figure  3.2:  Plot  of  seismic  moment  versus  The  numbers  refer 

to  the  event  list  of  Tables  I and  III.  The  standard  deviations 

on  the  M values  are  also  given, 
s 

Figure  3.3:  Plot  c.i  spectral  level  at  1 Hz  versus  ra^  for  two  assumptions 

of  the  shape  of  the  source  pulse  for  body  waves—rectangular 
(open  triangles)  and  Isosceles  triangular  (closed  circles). 
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4.  Computation  of  the  Source  Spectrum  from  a Crustal  Earthquake 

Gary  Lundquist 

Introduction 

We  showed  In  th’  last  semi-annual  technical  report  that  estimation  of 
the  power  spectrum  by  taking  the  Fourier  Transform  of  the  autocorrelation 
was  less  noisy  and  therefore  a more  desirable  estimation  than  direct 
application  of  the  FFT  algorithm.  We  continue  that  analysis  here  with  the 
objective  of  formulating  a practical  algorithm  for  computing  seismic  source 
spectra. 

This  study  is  divided  into  two  parts.  In  the  first  part,  we  extend 
the  autocorrelation  technique  for  power  spectral  estimation  to  the  estimation 
of  seismogram  amplitude  spectra.  In  particular,  we  examine  windowing 
techniques  which  sraooM'  the  spectrum. 

In  the  secend  part,  we  review  geophysical  corrections  to  the  seismogram 
spectrum.  In  part,  the  objective  is  to  understand  the  inherent  difficulties 
in  dealing  with  crustal  events  for  which  surface  reflections  and  multiple 
crustal  reverberations  are  part  of  the  P or  S wave.  We  also  examine 
anelastic  attenuation,  because  the  final  spectrum  depends  critically  upon 
the  model  used. 

Practical  Application  of  the  Autocorrelation  Technique  to  Estimation  of 
Seismogram  Spectra 

Power  spectral  estimation  from  the  autospectrum  requires  specifica- 
tion of  the  following  values 
1.  Total  data  length 

a.  Choice  of  first  and  last  points 

b.  Choice  of  sample  interval 

c.  Method  of  truncatii.g  the  raw  data 


A 

“1 


f 
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2.  Total  autocorrelation  length 

a.  Number  of  autocorrelation  lags  to  be  used  in  transform 

b.  Method  of  truncating  the  autocorrelation 

4 

We  will  test  various  choices  of  these  parameters  on  the  filtered  6 

k 

sequence  and  on  seismograms  for  the  Tien  Shan  earthquake  of  June  5,  1970. 
The  tests  will  be  conducted  in  the  context  of  the  general  flow  diagram 
for  autospectral  estimation  given  in  Figure  4.1. 

Data  Length  Testa 

We  will  commence  the  estimation  of  autospectra  with  raw  data  con- 
sisting ol  iionuniformly  sampled  scisoiogiaoi.';.  The  specification  ol  first 
and  last  points  must  be  baaed  upon  the  geophysics,  and  will  not  be 

discussed  here.  Yet,  given  a set  of  points  to  be  analyzed,  we  roust 
consider  how  we  remove  thor.c  points  from  the  rest  of  the  data  sequriuo. 

A simple  test  verifies  that  truncating  the  data  sequence  off  the 
zero  line  introduces  high  frequency  noise  into  the  autospectrun,  and  a 
sharper  truncation  creates  a greater  increase  in  high  frequency.  Since 
we  enn  decrease  the  high  frequency  content  1>y  increasing  rectord  length, 
this  behavior  is  certainly  related  to  truncation  and  not  to  data  length. 

We  cannot  apply  a smooth  data  window  to  the  raw  data  without 
seriously  distorting,  for  instance,  an  impulsive  P wave;  so  we  use  the 
following  method.  Given  a start  time,  we  ^earch  back  in  time  for  the 
first  zero  crossing  and  interpolate  a point  at  zero  amplitude.  That 
point,  then,  is  the  first  point  of  the  sequence  to  be  transformed,  and 
we  may  interpolate  a uniformly  sampled  sequence  considering  the  zero 
crocsing  to  he  at  time  zero.  The  end  of  ihr  trunciiind  sequence  may  he 
specified  by  desired  time  length  or  by  end  time.  In  either  case,  we  can 
apply  a taper  to  the  lest  few  points  rather  than  cicarrhing  for  anotiu'r 
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zero  crossing.  Testing  shows  that  a Gaussian  taper  applied  to  the  last 
ten  points  jf  the  interpolated  sequence  effectively  eliminates  the  effect 
of  truncating  off  the  zero  line.  To  summarize,  we  apply  a rectangular 
data  window  from  a point  interpolated  at  a zero  crossing  to  a point  a 
specified  length  away,  then  we  taper  the  end  of  the  data  smoothly  to 
zero,  leaving  the  data  virtually  undistorted. 

We  specify  the  sample  interval.  At  , according  to  the  frequencies 
we  want  to  resolve.  Sample  theorem  gives 

f -J-- 

Nyqulst  2At 


i 


1 

i 


where  the  Nyquist  frequency  is  the  highest  frequency  of  interest.  For 
the  case  of  seismograms,  the  choice  of  frequencies  has  been  made  for  us 
in  the  design  of  the  seismograph  systems.  The  systems  are  bandpass 
filters  with  half  power  points  (0.707  amplitude)  at  1.1  and  2.5  Hz  for 
WWSSN  short  period  and  0.035  to  O.IS  Hz  for  long  period  (15  - 100  Bystero). 
A more  conservative  bandwidth  for  the  choice  of  sample  Interval  is  the 
0.1  amplitude  point,  0.5  to  7.5  Hz  for  short  period  and  0.009  to  1.0  Hz 
for  long  period.  Thus  sample  intervals  of  0.067  and  0.5  sec  are  adequate. 

We  note  that  the  only  result  of  oversampling  is  the  extension  of 
the  Nyquist  frequency.  We  get^no  more  detail  in  the  pass  band  of  the 
instrument,  but  halving  the  sample  interval  doubles  the  number  of  fre- 
quencies in  the  spectrum,  and  all  of  the  extra  points  sample  frequencies 
wliii’li  nrc  .ibove  the  oM  NyqniHt  frequency.  Thus  ovrrnaniplin)’.  cokI;: 
money  to  produce  points  in  frequency  which  we  would  ignore  anyway  for 
being  well  outside  the  p.iKshand  of  the  instrument. 


i 


1 

i 

1 

1 
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Lag  Window  Comparisons 

Standard  methods  of  spectral  esti’:Atlon  from  the  autocorrelation  use 
a window  to  limit  the  number  of  autocorrelation  lags  used  in  the  Fourier 

I 

transform.  Given  a specific  autocorrelation,  we  have  two  parameters  to 
vary:  the  window  length  and  the  shape  of  l ha  window.  Of  course  tlu; 
autospectrum  and  the  perlodogram  are  Identical  If  we  vez  the  maximum 
autocorrelation  lag  length;  so  we  must  show  that  the  spectrum  Is  Improved 
by  windowing. 

We  will  consider  five  lag  windows  defined  as  follows: 
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The  rectangular  and  Gaussian  windows  are  endfonus  oi  a series  of 
windows  obtained  by  convolving  the  rectangular  window  with  itself.  The 
Bartlett  or  triangular  window  results  from  one  convolution,  nnd  a Parzen 
window  can  be  defined  as  the  result  of  four  convolutions.  As  the  number 
of  convolutions  increases,  the  central  limit  theorem  requires  the  resulting 
curves  to  approach  the  Normal  or  Gaussian  curve. 

The  transform  of  the  rectangular  window  is  a function  of  the  form 
sin  x/x  , soBStinas  called  a sine  function  or  the  Dirichlet  kernel.  The 
transform  of  the  Hanning  window  is  a linear  sum  of  shifted  Dirichlet 
keiii'.  ls  bocati-^c  ion  bv  rosinc  Tunc  t . ion  i:.  equivalent  l.o  a 

frequency  shift. 

M + 

where  D (u)  is  the  Dirichlet  kernel  of  an  N L ■*  1 poire  sequence. 

N ,, 

The  Hlackman  window  adds  a term  in  cos  , and  a generalized  Dirichlet 

la 

combination  nay  be  formed  as 

m 

(k)  - Aj^  co8^~  , |k|  1 L 

n*0 

Figure  ^ .2  compares  the  various  window  shapes  tested  for  this  st.udy. 

The  truncation  of  the  autocorrelation  in  time  by  a window  is 
equivalent  to  convolution  in  the  frequency  domain  of  the  complete  auto- 
spectrum with  the  transform  of  the  window  function.  A desirable  window 
is  one  which  Irunc.ilcs  the  autocorrelation  function  5:moothly  and  provides 
adequate  noise  reduction  in  the  autospectrum. 
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The  rectangular  window  has  a frequency  donain  character  which 
depends  upon  the  window  length  relative  to  a power  of  two.  To  an  FFT 
algorithm  a rectangular  window  of  lengthy  equal  to  a power  of  two  looks 
like  an  infinite  sequence  of  I's,  and  the  transfore  is  correspondingly  a 
spike  at  zero  frequency.  The  effect  is  no  smoothing  and  a noisy  auto- 
spectrun.  If  the  length  is  one  less  then  a power  of  two,  the  window 
appears  to  be  the  sum  of  a constant  and  a spike  of  amplitude  minus  one. 

The  result  in  frequency  is  the  sum  of  a spike  at  D.C.  and  a constant 
amplitude  spectrum.  For  any  other  length,  the  rectangular  window  trans- 
forms to  a sine  function  as  predicted  by  continuous  transform  theory. 
Depending  upon  the  window  length  and  the  smoothness  of  the  true  spectrum, 
the  odd  length,  one  less  than  a power  of  two,  is  likely  to  give  a 
rectangular  window  frequency  response  much  better  than  other  lengths 
which  give  a sine  function  response. 

The  Hanning  window  also  transforms  differently  for  different  lengths. 
In  particular,  the  Hanning  window  is  optimized  if  the  lag  length  is  one 
less  than  a pswer  of  two.  For  that  case,  only  the  firat  two  frequency 

domain  points  have  any  power,  ar  all  other  samples  have  amplitudes  below 
“26 

10  . Considering  the  sysmetry  of  toe  Fourier  transform,  we  sea  that 

this  means  the  Hanning  smoothing  is  a 3 point  weighted  running  average. 

In  terms  of  application,  this  suggastu  we  compute  the  rutocorrelatlon  to 
lag  length  equal  to  a power  of  two.  Hien  upon  folding  (sec  Figure  4*3), 
we  have  a length  one  less  than  the  next  power  of  2 and  we  get  the 
optimum  Hanning  ''csponac. 


t 


89 


The  Hanming  window  is  very  similar  to  the  Hanning  except  for  the 
special  length  cases.  In  general.  It  shows  higher  resolution  but  a 
slower  decay  rate.  For  the  length  one  less  than  a power  of  two,  the 
Hannning  response  is  similar  to  that  of  a reutangular  window  in  that  it 
att.*yins  a finite  value  and  remains  constant. 

The  response  for  the  Bartlett  window  has  been  examined  and  does  not 
have  any  particular  advantages. 

The  Gaussian  window  differs  from  the  Hanning  type  In  that  It  has 
no  reference  to  window  length  In  Its  formula  and  that  It  depends  on  a 
parameter  u'  whicli  relates  to  the  width  o..  the  window.  The  Hanning 
window  for  256  points  Is  just  a scaled  up  version  of  the  64  point  Hanning 
window.  The  Gaussian  window,  on  the  other  hand,  uses  exactly  the  same 
amplitudes  In  each  window  (for  constant  m)  as  far  as  the  lengths  allow. 
Thus,  by  going  96  points  beyond  the  32  point  cutoff  for  a 64  point 
window,  the  256  point  Gaussian  window  (w  «=  0.07)  attains  amplitudes  as 
low  as  0.00005  and  smooths  the  last  40  points  quite  severely.  Obviously 
the  "width"  parameter,  u , must  be  varied  with  the  lag  length  to  get 
equivalent  responses  for  different  window  lengths.  For  u very  small 
(0.067  for  L = 64)  the  Gaussian  window  approximates  u rectangular  window, 
but  as  a)  Increases,  the  window  becomes  narrower  and  more  bell  shaped. 

The  next  step  in  evaluating  the  windows  Is  to  apply  them  to  a signal 
with  a known  spectrum,  In  this  case  the  frequency  response  of  a Butterworth 
lowpass  filter.  The  theoretical  spectrum  (perlodogram)  is  characterized 
by  a continuous  and  rapid  rolloff  beyond  the  cutoff  frequency  at  0.2  Hz. 
Each  of  the  autospectra,  however,  sliows  n lower  limit  lo  the  decay.  In 
this  context,  a good  window  is  one  which  yields  an  outospectnim  which  if  a 
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close  approximation  to  the  theoretical  spectrum.  Such  a window  minimizes 
the  amplitude  of  the  noise  which  Interrupts  the  high  frequency  rolloff 
and  one  which  gives  the  theoretical  decay  rate  down  to  Its  noise  limit. 

I 

Both  the  periodogram  and  the  Dfrlchlet  smoothed  autospsctrum  decay 
to  10  ^ at  0.3  Hz,  so  the  second  criterion  Is  satisfied.  However,  the 
decay  falls  to  continue  after  0.32  Hz,  and  the  final  level  at  0.5  Hz  Is 
3 X 10  Applying  a Hanning  window  maintains  the  character  of  the 
rolloff  to  0.4  Hz,  and  gives  a final  level  of  8 x 10  or  2 1/2  ordsrs 
of  magnitude  better.  The  Hamming  giv'.j  a higher  noise  level  than  the 
haniiing,  and  is  only  slightly  better  than  the  retLangular  window  for 
window  length  equal  to  a power  of  two. 

For  a fixed  window  length,  the  Gaussian  smoothed  autospectra  depend 
on  (ii  . For  small  lo  rhe  response  is  that  given  by  a rectangular 
window,  with  correct  rolloff  but  high  noise  level.  As  u Increases,  rhe 

noise  level  decreases,  but  the  rolloff  is  distorted.  Thus,  at  id  = 0.32, 
the  noise  level  is  lower  than  that  for  the  Hanning  window,  but  the  level 
10  ^ Is  not  reached  until  after  0.32  Hz,  and  for  id  larger  than  0.2, 
the  distortion  of  the  slope  is  unaceeptably  poor. 

The  Bartlett  window  has  an  unacceptably  high  noise  level. 

The  observations  made  so  far  apply  to  an  extremely  simple  ideal 
spectrum,  and  we  found  the  Hanning  window  to  be  the  best.  We  now  form 
a window  comparison  test  for  earthquake  seismograms.  Figure  4.3  shows 
the  smoothed  levels  of  the  mitospeetra  relative  to  the  periodogram  for 
BAG  LPZ  and  BAG  SPZ  for  the  event  of  June  5,  1970.  The  conclusion  for 
the  Long  period  Lost  is  obvious,  iiic  Hanning  window  h.is  the  lowest  noise 
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level  and  Its  slope  follows  the  pv'rlodogram  nicely.  Also,  though  net 
shown  In  this  figure,  the  Hanning  s.noothed  autospectrum  has  the  definite 
appearance  of  a spectrum,  while  the  rectangular  and  Hamming  spectra  are 
dominated  by  window  noise.  The  Gaussian  spectra  (m  “ 0.101)  Is  smoother 
than  the  Hanning,  but  It  does  not  have  the  correct  slope  In  the  range 
0.15  - 0.3  Hz. 

The  short  period  comparison  Is  not  so  easily  analyzed.  The  Hanning, 
Hammier.,  Gaussian  (ca  = 0.081)  and  Bartlett  windows  all  give  nearly 
equivalent  responses.  The  Hanning  spectrum  has  a more  nearly  correct 
decay  rule,  while  the  Hamming  has  slightly  more  resolution.  Since  the 
long  period  data  so  clearly  favors  the  Hanning  window,  we  will  use  the 
Hanning  window  In  all  further  tests. 

Window  Closing  Tests 

We  can  define  the  bandwidth  of  a lag  window  In  several  ways.  In 
terms  of  filter  theory,  the  bandwidth  Is  the  distance  between  half 
power  points.  In  statistical  applications,  the  bandwidth  may  be  defined 
as  the  width  of  a rectangular  window  with  the  seme  variance^  But  whatever 
definition  we  use,  there  is  an  inverse  relationship  between  window  length 
in  time  and  bandwidth  in  frequency.  In  a window  closing  test,  we  compute 
smoothed  spectral  estimates  for  increasingly  wide  time-domain  windows, 
effectively  closing  the  frequency-demain  window.  If  we  visualize  the 
convolut.ion  process  in  frequency  as  a weighted  running  average,  then  it 
is  obvious  that  a narrower  bandwidth  will  resolve  more  of  the  detail  in 
the  original  spectrum.  Thus,  as  we  increase  the  number  of  lags  used,  we 
Improve  the  resolution  in  frequency  untii.,  in  the  limit,  we  obtain  the 
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perlodogram.  Lauklng  at  It  anotkei  way,  as  the  bandwidth  decreases,  the 
spectral  estimate  Is  a polynomial  of  Increasingly  high  degree  In  cos  2Trf, 
so  It  Is  easy  to  produce  spurious  peaks. 

The  question  we  must  answer  here  Is  when  to  stop  the  process  of 
narrowing  the  bandwidth.  The  answer  will  depend  upon  experience,  the 
amount  of  detail  desired,  and  personal  preference.  Sometimes  the  spectrum 
may  be  stable.  In  which  case  all  of  the  detail  will  appear  and  decreasing 
bandwidth  will  resolve  nothing  more.  But  more  often  the  spectmm  will 
become  unstable  and  noise  will  begin  to  mask  significant  detail.  The 
tlicorct  ii  .'H 1 filter  r«  onsc  (a  filtered  5^^  sequence)  is  a member  of  the 
first  class,  and  seismograms  belong  to  the  second  class. 

In  Figure  wc  show  the  window  closing  experiment  for  our  BAG  SPZ 

seisraogT:»ni  in  terns  of  amplitude  spectra.  The  number  of  lags  varies  from 
52  of  the  sequence  length  at  the  bottom  to  252  at  the.  top.  The  Fourier 
amplitude  spectrum  is  shown  at  the  top  for  comparison.  An  overlay  would 
show  that  the  slope  and  levels  are  defined  for  the  52  case,  and  further 
lag  length  only  adds  detail  to  the  spectrum.  As  mentioned,  the  amount 
of  detail  (or  noise)  which  Is  desliable  Is  a matter  of  personal  preference. 

f 

We  find  202  of  the  sequence  length  to  be  a maximum  and  102  to  be  preferable. 
For  shorter  time  sequences,  though,  a longer  transform  length  contributes 
to  the  smoothness  of  the  spectrum,  so  we  will  compromise.  For  short 
period  seismograms  we  wll*  use  152  of  the  complete  autocorrelation 
length  In  computation  of  t.he  autospectrum. 

Figure  A. 5 shows  the  amplitude  autospectra  of  the  window  closing 
experiment  on  BAG  l-l’Z.  Once  again,  a 20Z  lag  length  provides  the  basic 
spectral  shape  for  the  frequencies  above  1 Hz,  but  for  the  frequencies 
below  1 Hz,  wo  must  consider  the  autospectrum  technique  In  more  detail. 
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The  autocorrelation  function  has  the  following  properties:  it  is 

real  and  has  even  symmetry;  its  maxlmuni  occurs  at  zero  lag;  it  contains 
no  information  on  phase  relationships;  and  it  contains  the  same  ampli- 
tude information  as  the  original  signal.  If  a perludlclty  is  present 
in  the  time  sequence,  it  is  present  with  the  same  amplitude  but  zero 
phase  in  the  autocorrelation.  This  last  property  implies  that  resolution 
of  long  periods  requires  a long  length  of  autocorrelation.  Resolution 
of  a 30  sec  periodicity  requires  30  see  of  lag. 

Going  back  to  the  long  period  window  closing  test,  we  note  that  the 
/O,,  lag  length  resolves  a maximum  pciiodicity  of  1.9.?.  sec,  or  0.05?  Hjj, 
for  a data  length  of  320  points,  while  the  40Z  lag  length  resolves 
47  see.  Thus  the  apparent  Improvement  in  resolution  of  the  peaked 
character  is  real,  and  very  important.  We  must  compromise  between  the 
smoothness  of  the  spectnua  and  resolution  of  low  frequencies.  In  our 
opinion,  no  slgniflctnt  detail  is  resolved  for  lag  lengths  greater  than 
40%,  and  the  noise  level  is  still  acceptable  at  40Z  for  this  particular 
earthquake.  However,  an  unavoidable  conclusion  is  that  the  best  seismic 
spectra  will  require  a window  closing  experiment  for  the  long  period 
autospectra  of  each  individual  earthquake. 

Discussion 

We  have  examined  the  practical  problems  of  estimating  seismic  source 
spectra  up  to  the  point  of  applying  instrument  response  and  attenuation 
corrections.  We  have  chosen  a rectangular  data  window  with  a short 
taper  at  the  end,  and  we  found  the  Hanning  window  to  be  the  best  of  the 
lag  windows.  The  long  and  short  period  sample  intervals  should  bo  Icr.r; 
than  0.5  and  0.06  sec,  respectively,  and  the  lag  length  should  be  near 
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40Z  of  the  data  length  for  long  period  records  and  IbX  for  short  period 
records.  Finally,  although  we  cannot  specify  the  data  lengths,  experience 
has  shown  us  that  we  need  short-period  data  lengths  of  15  to  60  sec  and 
long-period  data  lengths  of  40  to  180  sec. 

The  flow  diagram  in  Figure  4.1  gives  the  steps  required  to  obtain 
power  autospcctra.  The  amplitude  aiitospectra  are  obtained  by  simply 
taking  the  square  root  of  each  frequency  sample  multiplied  by  the  number 
of  samples  in  the  autocorrelation,  then  multiplying  by  the  sample  Interval. 

The  results  of  our  tests  demonstrate  that  we  can  obtain  less  noisy 
and  thus  inoie  reliably  interpreted  seismic  spectra  using  the  autocorrelation 


technique. 
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Ceophystcal  Corrections  to  Seismogram  Spectra 

The  computation  of  a seismic  source  spectrum  from  the  estimates  of 

seismogram  spectra  at  several  stations  requires  that  a number  of  corrections 

be  applied.  We  shall  refer  to  a spectrum  which  as  been  corrected  for 

instrument  response  as  a station  spectrum.  The  corrections  which  Liust  be 

applied  to  station  spectra  to  obtain  a source  spectrum  are  intended  to 

invert  the  transmission  filter  in  the  Earth.  As  a first  approximation, 

we  may  write  the  station  spectrum  of  a body  wave  as  the  product  In  frequency 

of  Independent  effects  of  the  source  spectrim  F(ai)  with  transfer  functions 

of  the  source  Jirectivity  function,  oi  radiation  pattern,  R(0)  , nourco 

crustal  response,  C (i  , u)  , geometric  spreading,  G(A)  , anelastic 

s o 

attenjiation,  Q(A,  u)  , and  receiver  crustal  response,  Co(i,  , to)  . Here 

R n 

A is  distance,  measured  in  degrees  ot  central  angle  subtended  by  the  ray; 

0 is  azimuth,  i is  take-off  angle  from  the  source.,  and  1.  is  incident 
o n 

angle  at  the  base  of  the  receiver  crust.  That  is 

D(w,  A,  0)  = F(<o)  • R(0)  • Cp(ip,  u)  • G(A)  • Q(A,  to)  • eo)  (Bl) 

where  D is  the  station  spectrum  or  ground  displacement  spectrum.  Since 
the  object  is  the  source  spectrum  F , w;e  must  consider  correcting  for  each  of 
the  contributions  to  the  Earth  filter. 

Source  Crust 

We  wULL  flt.1  the  evaluation  of  the  source  crustal  response  to  be 
simpler  if  we  consider  constructing  a synthetic  seismogram  from  the  multiple 
contributions  of  waves  reflected  internally  in  the  crust,  rather  than 
trying  to  decompose  an  observed  seismogram.  A critical  point  in  this 
entire  study  is  that  we  are  dealing  with  crustal  events  for  which  the 
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Observed  P vrave  is  not  Isolated  from  surface  reflections  or  multiple 
reverberations  in  the  source  crust. 

In  a separate  part  of  this  study,  Herrmann  (1975)  succeeded  in  syn- 

* * 

theslzlng  long-period  P-wave  seismograms  of  crustal  events.  To  do  so, 
he  used  the  four  basic  waves  which  can  contribute  to  the  teleselsmlc  P 
wave  and  added  31  adJltlonal  waves  from  .cxnistal  reverberations  near  the 
source.  In  Figure  4.6  we  show  the  shapen  and  relative  amplitudes  of  the 
contributions  of  those  four  waves  and  their  individual  reverberatlrns  to 
the  synthetic  seismogram  of  a purely  strike-slip  cource  at  a station  with 
11,  l.'if  ivu  ar.fiuuth  of  40''  and  distance,  of  45*’.  The  dovrugoing  or  ddifcr 
P wave  la  reflected  at  the  crust-mantle  boundary,  producl’.ig  upgolng  P 
and  SV  waves.  Kach  of  those,  arc  reflected  at  the  surface  with  twofold 
convf.rsfou,  and  so  on~  The  signal  pictured  In  Figure  4.6  is  the  contribution 
at  the  receiver  of  this  sequence  of  waves.  Similarly,  the  contributions 
from  upgoing  p wave  (pP)  and  up  and  downgoing  SV  (sP  and  S„P)  are  presented. 

The.  obvious  feature  is  the  much  larger  amplitude  of  sP  relative  to  P. 

T.f  the  source  radiated  equal  P and  S’v  energy,  theoretical  reflection 
and  transmission  coefficients  would  give  a P amplitude  at  thl.s  particular 
receiver  more  than  three  times  as  large  as  the  sP  amplitude.  The  observed 
amplitude  ratio  results  from  the  efficiency  of  a shear  dislocation  in 
radiating  shear  energy.  Of  course  the  observed  amplitude  ratio  will  vary 
from  zero  to  Infinity  depending  upon  station  location  and  focal  mechanism. 

But  over  most  of  the  focal  sphere  the  contribution  '"o  the  P^ave  seismogram 
of  a crustal  event  from  the  upgolng  SV  will  be  greater  than  or  equal  to 
the  contrtbullon  from  the  direct  P.  Tn  other  words,  if  we  Ignore  the 
interference  phenomenon  for  a momeitt,  the  i-esponse  of  the  source  crust  is 
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to  In.crea'if.  the  apparent  spectral  level,  by  a factor  depending  upon  focal 
mechanism,  over  the  level  we  would  measure  for  an  Isolated  P-wave. 

The  question  we  must  answer  Is  whether  we  can  correct  fo:  the  source 
crustal  response.  This  leads  us  to  consider  the  Interference  of  the  four 
wave  types.  Interfercuice  depends  upon  relative  time  delays  which  are 
functions  of  crustal  thickness  and  source  depths,  rather  than  focal 
mechanisms . 

The  contribution  of  each  wave  type  may  be  written  as  a Sian  of  roflected 
and  refracted  components 

K 

f(t)  - 2 f^(t  - T^)  (1) 

1-1 


where  Is  the  source-time  function  for  the  particular  wave  type 

being  modelled;  the  are  amplitudes  determined  by  reflection  and 

refraction  coefficients  and  by  the  radiation  pattern,  and  the  are 
delay  time  associated  with  the  individual  components.  We  can  sum  over 
the  four  wave  types  and  compute  the  spectrum  of  the  total.  That  Is,  if 
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wheio  f^(t)  is  Lhc  source-time  funct.ioo  for  douiij-olag  P waves; 
corresponds  to  upgolng  P waves,  and  and  to  upgoing  and  down- 

going S waves,  respectively.  We  may  simplify  (2)  by  making  assumptions 
about  these  source  time  functions.  If  we  neglect  rupture  effects,  then 
we  may  assiane  f^(t)  = f^^(t)  - fp(t)  , and  f^(t)  -=  f^(t)  - f^(t)  . Further, 

if  we  model  the  source  as  a iouble  couple  point  source,  then  f (t)  “ 

P 

f (t)  “ f(t)  . In  that  case 
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where  the  range  of  n includes  all  of  the  possible  contributions. 

We  can  see  from  (3)  that  the  lopo-period  (u=«o)  level  of  tt.e  total 
spectrum  is  just  the  sum  of  the  levels  of  the  spectra  of  the  Individual 
contributions.  Thus  the  computation  of  seismic  moment  is  seriously 
affected  if  we  can  actually  estimate  the  DC  term  in  the  spectrum. 

For  any  other  frequency,  the  spectral  level  depends  upon  che  inter- 
ference of  the  sinusoidal  modalation  caused  by  the  delay.  Thus  we  cannot 
remove  the  effect  of  the  source  crust  in  any  simple  manner.  The  only 
technique  which  seems  possible  is  Herrmann's  method  of  modelling  the  P 
wave  in  the  time  domain  as  a sura  of  waves.  Once  a fit  has  been  achieved, 
then  one  possible  form  of  the  contribution  of  any  individual  wave  is 


known. 


It  Is  appropriate  to  note  here  tluL  we  ahould  not  arbitrarily  nnaunu*. 
that  P and  S wive  source-tine  functions  are  the  same.  If  they  differ, 
then  the  total  DC  level  still  varies  as  a sura  of  DC  levels.  However, 
if  P and  S comer  frequencies  differ,  then  a corner  frequency  of  the  total 
spectrum  should  reflect  that  difference,  for  a crustti  event,  by  being 
iQQi>g  rounded  or  perhaps  two  cornered,  rather  than  being  a single  sharp 
band. 

To  suamarlze  briefly,  we  should  expect  the  spectrum  of  a P wave 
generated  by  a crustal  source  to  have  a higher  long  period  level  than 
that  y’ciicrated  by  the  same  focal  mecbiinipm  at  depth.  We  should  expect 
such  a spectrum  to  show  the  interference  of  all  of  the  crustal  rever- 
berations and  in  particular,  of  the  direct  P and  the  sP.  This  effect  is 
to  apply  a sinusoidal  modulation  to  the  sum  of  the  separate  spectra 
and  is  an  effect  not  simply  removed.  And  finally,  if  the  source  time 
functions  of  the  various  contr.^butlons  are  all  identical,  we  should 
expect  any  cornero  in  the  spectrum  to  be  sharp  and  well  resolved.  If 
such  corners  are  rounded  or  poorly  resolved,  it  may  be  evidence  that  the 
source  does  not  appear  exactly  the  same  to  an  S wave  es  it  does  to  a P 
wave. 

Radiation  Pattern 

In  the  last  section  we  showed  that  the  contribution  from  the  four 
wave  types  could  not  be  separated  by  decomposition  of  observed  tele- 
seismic  ground  motion  due  to  a crustal  earthquake.  Thus  the  question  of 
application  of  a radiation  pattern  correction  to  such  a record  reduces 
tc  the  question  of  whether  or  not  we  y apply  one  correction  for  all 
contributions  to  the  P wave.  The  answer  is  obviously  no.  In  general 
the  downgoing  P and  the  upgoing  SV  will  sample  different  radiation 
patterns,  and  we  cannot  expect  to  correct  the  entire  telcseismic  P wave 
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by  radiation  pattern  computed  for  P alone. 

We  can  partially  correct  for  radiation  pattern  by  averagings  but 
even  a perfect  average  of  a radiation  pattern  of  naxlmua  amplitude  of 
one  la  not  equal  to  one.  If  we  assume  a perfect  average,  then  we  should 
correct  the  P wave  by  1/0.424  and  the  S wave  by  1/0.593. 

Attenuation 

We  can  separate  amplitude  loss  along  the  transmission  path  into 
a frequency  dependant  part,  the  anelaatlc  attenuation,  and  a frequency 
Independent  part,  the  geometric  attenuation,  or  spreading.  In  the 
computation  of  earthquake  nagnltudci;  ■ wc  use  Gutenberg  and  Rich'rer's  (1956) 
Q to  correct  for  both  types  of  attenuation,  aasunlng  no  frequancy  de- 
pendence. For  the  correction  of  spectra,  however,  we  must  apply  separate 
factors. 

We  take  our  geometric  attenuation  correction  from  Duda  (1971). 

He  notes  that  a basic  falling  of  strictly  empirical  attenuation  factors 
Is  the  inability  to  relate  amplitudes  of  earthquakes  from  different  depths. 
He  '.hen  evaluates  Bullen's  (1953)  equation  for  amplitude  as  a function  of 
distance  in  a spherically  synmetrlc  Farth.  The  travel  path  Is  defined  by 
ray  theory  and  by  the  velocity  structure  of  Herrin  et  al.  (1968).  In  this 
model,  no  velocity  Inversions  are  allowed,  and  the  resulting  corrections 
for  distance  and  depth  arc  much  smoother  than  Richter's.  Nuttll  (1972) 
computed  geometric  attenuation  empirically,  using  nuclear  explosion 
sources.  In  his  comparisons,  Nuttll  found  his  and  Duda's  curves  to  be 
very  similar  in  the  distance  range  25*  — 95*,  and  concluded  both  curves 
were  preferable  tc  Richter's  Q for  magnitude  determination. 

The  model  for  anelastic  attenuation  is  more  Involved.  In  a linear, 
perfectly  elastic  medium,  the  amplitude  of  a propagating  stress  wave  is 
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where  A is  amplitude  and  r is  distance  from  the  source.  If  the 
medium  is  linear,  but  viscoelastic,  then  the  wave  number,  k , may  be 
considered  complej;.  Thus,  while  there  is  no  absorbtion  in  an  elastic, 
medium,  the  amplitude  in  a viscoelastic  medium  damps  according  to 


Ae 


-yr  + l(kr  - mt) 


where  k is  the  real  part  and  y the  imaginary  part,  respectively 

of  the  complex  wave  number,  y is  called  the  spatial  attenuation  coefficient 

md  is  defined  by  (Nutt  Li,  1973) 
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where  Q is  the  dimension"'  ss  quality  factor  which  describes  the 
fractional  elastic  energy  dissipated  per  cycle,  V is  the  wave  velocity, 
f is  frequency  and  P is  period.  A region  of  high  Q is  a region  of 
low  attenuation,  and  a low  Q region  absorbs  relatively  more  energy 
from  the  passing  wave. 

Since  Q is  a per  cycle  charactlzation  of  absorbtion,  it  is  con- 
venient to  integrate  Q along  the  ray  path,  C . 
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Then  the  total 


Q-type  atenuatlon  along  the  ray  path  is  given  by 


-to  / 


2Q 


(7) 


The  numerical  integration  requires  a velocity  model  to  obtain  the  travel 
time,  T , and  a Q model  to  obtain  the  effective  Q , Q ff  • We 
have  also  assumed  Q to  be  Independent  of  frequency.  That  is,  the 
integration  Is  performed  by  ;:indi£ig  the  number  of  cycles  (ur.ing 
travel— time  and  ’.iiocity),  and  then  summing  the  same  Q for  each  cycle, 
regardless  of  wavelength  or  frequency. 

Currently  usee  Q models  show  a ^one  of  high  Q , the  lithosphere, 
ever  lying  a region  of  lov#  Q , the  asthenosphere,  which  in  turn  overlies 
the  re  st  of  the  mantle  which  is  also  a high  Q region.  MoeL  of  the 
attenuation  of  seismic,  waves  appears  to  take  place  in  the  asthenasphere. 
The  result  of  the  absorption  being  concentrated  in  a thin  zone  is  that 
t-otal  absorption  ts  a function  of  distance  is  nearly  constant  because 
for  all  teleseismic  distances,  all  rays  spend  nearly  equal  times  in  the. 
asthenosphere.  The  mantle  below  the -vsthenosphere  is  approximately  non- 
absorbing, with  Q large  and  Increasing  with  depth.  This  accounts  for 
the  observation  that  the  ratio  Is  u«.arly  constant  with  re.spcct 

to  epicentral  distance,  though  travel  time  and  effective  Q Increase. 

The  commonly  observed  ratio  for  P waves  is  - 1.0  , and  Anderson 

et  al.  (1965)  show  that  the  ratio  of  Q for  S waves  to  Q for  P waves 
uhould  be  2.25,  by  requiring  density  to  be  real  and  by  Ignoring  losses 
due  to  Imperfect  inertia. 

Solomon  (1971)  has  reviewed  and  supplemented  a large  body  of  evidence 
for  lateral  variation  of  Q , and  moreover,  he  presents  evidence  in 
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favor  of  a Q model  which  depends  upon  frequency.  This  is  in  agreement 
with  Tsai  and  Aki  (1969)  who  found  a Q model  derived  from  surface  waves 
to  be  unacceptable  for  correcting  hl^^h  frequency  (f  > .1  Hz)  waves. 

The  mechanism  proposed  by  Solomon  for  anelastlc  attenuation  is 
viscous  damping  in  thin  films  of  fluid  in  a partially  molten  asthenosphere. 
We  can  visualize  such  a medium  as  a concentration  of  dislocations,  each 
capable  of  moving  under  an  applied  stress.  In  such  a medium,  the  shear 
modulus  is  also  a function  of  frequency.  Vary  little  energy  is  required 
to  force  long  periot’  relative  motions  of  particles  when  the  grains  are 
lubricated  by  a partial  melt.  That  is,  the  shear  modulus  is  very  low  and 
little  energy  is  absorbed  by  the  medium.  At  very  high  frequencies,  the 
viscosity  is  such  that  the  medium  absorbs  very  little  energy  because  the 
effective  stiffness  is  very  high.  At  some  intermediate  frequency,  the 
energy  absor'ljtion  is  a maximum,  and  we  can  refer  to  the  behavior  in  the 
neighborhood  of  this  frequency  as  a relaxation.  That  is,  the  medium 
achieves  a new  equilibrium  under  an  applied  stress  in  a relaxarion  time, 

T , characteristic  of  the  medium.  When  x “ 1/u  , then  the  energy 
absorbed  by  the  medium  is  maximum,  and  Q is  a minimum.  Since  the  loss 
is  in  shear,  this  model  agrees  with  the  relatively  larger  attenuation  of 
S waves  versus  P waves. 

The  functional  dependence  of  Q on  frequency  caused  by  a particular 
relaxation  phenomenon  is  given  by 

1 Au  “>  T 

i » Vi"  . M2  T ---;  ,2 

Q (1  - Ap)  1 + 

where  Ay  is  the  change  in  the  shear  modulus.  A graph  of  this  function 
is  a bell  shaped  curve  with  1/Q  a maximum  at  ux  = 1 . Liu  and 
Archambeau  (1975)  found  a Solomon  type  model  to  be  the  best  of  those 
tested  against  observed  attenuation  of  free  oscillation  data. 
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We  can  reunrite  our  definition  of  y with  a frequency  dependent 

Q , 

^ 

The  prospect  of  computing  a separate  for  each  frequency  at  each 

distance  is  not  pleasant  and  we  propose  the  following  simplification 

Q(*.  f)  " qW 

JIici'  the  ill  t.entirtirii  cxpcnent  is 

y ' “Tn^-  R(f)  (10) 

In  otiier  words,  under  this  assumption  we  may  still  use  constant  T/0 

^ef  f 

ratios  and  simply  modulate  the  exponent  with  Solomon’s  relaxation 
fund  ion. 
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^^ef  f (1  - 1 + (o,t)^ 


(11) 


Here  the  dependence  on  Ap  becomes  a scale  factor  on  the  ratio  T/C 

^eff  ' 

Since  the  last  factor  has  a maximum  value  1/2  at  ut  » 1 , we  let  the 

scale  factor  equal  2.0.  Any  other  choice  of  scale  factor  is  equivalent 
to  changing  the  T/Q^^^  ratio. 

**e  will  attempt  to  choose  a Q function  on  the  basis  of  observed 
spectra.  In  this  choice,  we  require  only  that  the  spectrum  exhibit 
some  sort  of  comer  beyond  which  the  high  frequency  decay  approaches 

3 ~ 

1/f  dependence  asymptotically.  The  choice  of  1/f^  dependnence 


agree  with  m - m.  data  and  energy  finiteness  arguments.  By  specifying 
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that  we  see  a 1/f  dependence  within  the  passbands  of  our  Instruments, 
we  are  assuming  a moderately  high  to  high  rupture  velocity. 

Figure  4.7  shows  a few  representative  modulations  of  the  total 
frequency  dependent  anelastlc  attenuation  factor  as  a function  of  the 
relaxation  time,  t . Note  that  as  t Increases,  the  peak  of  the 
relaxation  function  moves  to  lower  frequencies,  causing  reduction  In 
the  amount  of  attenuation  Implied  for  higher  frequencies.  Indeed,  If 
we  plotted  the  curve  for  t - 2 sec  , which  LLi  and  Archambeau  used,  we 
would  03 serve  almost  no  attenuation  at  all.  Their  requirements  were  for 
modulation  at  .cry  low  frequencies,  and  thelt  data  had  no  resolution 
at  periods  less  than  20  sec. 

Figure  4.8  displays  the  effect  of  the  relaxation  modulation  on  the 
average  P-wave  spectrum  for  the  Tien  Shan  event  of  February  2,  1969. 

At  the  top  is  the  spectrum  corrected  with  an  unmodulated  T/Q  = ]. 0 . 
For  ( ■=  0.2  sec  , the  short  and  long  period  decay  rates  are  equal, 
giving  a well  defined  corner  at  .1  Hz.  For  0.2  < t the  peak  of  the 
relaxation  curve  occurs  at  frequencies  below  .8  hz,  and  the  spectrum  shows 
two  distinct  corners.  That  is,  if  a relaxation  phenomenon  with  t in 
this  range  controls  the  anelastlc  attenuation,  then  the  true  two-corner 
character  of  an  earthquake  spectrum  would  be  masked  by  the  transmis.sion 
A*’  "t  increases  beyond  0.25,  the  corners  on  the  spectrum 
slowly  converge.  At  the  same  time,  the  spectrum  approaches  the  station 
spectrum,  and  Indeed  the  station  spectrum  has  the  desired  1/f^  slope. 
However,  the  amount  of  total  attenuation  predicted  when  t is  larger 
than  0.5  is  unrealistically  small  at  high  frequencies.  For  the  purposes 
of  this  report,  we  will  chose  t = 0.3  as  the  appropriate  parameter 


for  the  relaxation  model. 
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Further  modelling  is  required  to  test  the  physical  significance  of 
this  choice.  We  have  chosen  a ratio  vhlch  is  constant  with 

distance,  'ihis  Implies,  of  course,  an  Increasing  Q with  depth. 

Thu  modulation,  however  probabl  y applies ‘only  to  the  asthenosphere 
because  of  the  physics  of  the  relaxation  model,  and  is  constant  with 
distance  only  because  all  teleselsmlc  body  wave  spend  nearly  the  saiD2  time 
In  the  asthenosphere. 

Solomon  tested  models  with  more  than  one  relaxation  peak,  and  the 
real  Farth  probably  la  a multiple  relaxation  medium.  Additional  ttsts 
wlJI.  be  run  tii  relate  the  physical  properties  of  the  relaxations  to  the 
expected  characteristics  of  out  seitimlc  source  spe.:tra. 

' Lver  Crust . 

lien-Menahem,  Smith  and  Tcng  (1965)  implemented  Haskell's  layer  matrix 
technique  for  six  varied  crustal  models.  Their  computed  vertical  P- 
wave  crustal  amplitude  responses  show  gains  of  2.0  - 2.6  for  all  but  one 
model  (incident  angles  about  30*  and  incident  waves  in  the  period  range 
1 sec  - 40  sec).  Over  a period  range  1 - 100  sec,  we  can  approximate 
the  receiver  crustal  response  by  a gain  of  2.3  with  a maximum  error  of 
.7  at  low  frequencies  which  tend  to  8«ie  the  layered  system  as  a half 
space . 

A refinement  which  should  be  considered  for  this  study  Is  application 
of  receiver  crustal  corrections  for  each  station.  The  result  of  Ben- 
Menaheu  ^ al.  show  that  the  various  crustal  models  tested  give  similar 
responses  for  periods  greater  than  10  sec  but  that  the  high  frequency 
responses  vary.  Since  such  corrections  are  not  available,  we  apply 
a stralght-llnc  correction  and  assume  the  average  spectrum  will  not  be 
distorted  by  individual  errors. 
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S'immary 

We  have  found  each  of  the  corrections  Implied  by  equation  B1  to  be 
complicated  functions  of  various  combinations  of  the  parameters  epicentral 
distance,  soiii''e  depth,  station  azimuth,  the  Earth's  velocity  and  Q 
structures,  and  frequency. 

We  detenained  that  the  source  crust  tends  to  increase  the  long  period 
level  of  the  spectrum  but  that  the  interference  of  various  waveforms  could 
not  brt  undone  for  the  spectrum  as  a whole.  The  source  crustal  response 
also  depends  upon  orientation  of  the  focal  mechanism,  and  our  inability 
to  correct  for  crustal  response  means  that  we  also  cannot  correct  for 
radiation  pattern  because  the  various  contributions  do  not  always 
come  from  similar  radiation  patterns.  Rather  tha/'.  compute  an  intricate 
and  perhaps  viong  conextlon  factor,  we  hope  to  mlnlTtize  the  maximum 
error  by  setting  C (i,  oj)  =1  ■ R(0)  = 1 and  averaging  over  several 

o 

stations.  The  average  over  a perfectly  sampled  radiation  pattern, 
however  gives  .A94  for  P or  .593  for  S . Thus  we  divide  the  average 
P or  S wave  spectra  by  the  appropriate  number  on  the  assumption  that  our 
average  value  is  lower  than  the  maximum  of  the  dominant  radiation  pattern. 

The  geometric  attenuation,  or  spreading  correction  given  by  Duda  was 
chosen.  A typical  value  for  epicentral  distance  of  45“  is  ('''5  x 10^’). 

The  anelastlc  attenuation  turns  out  to  be  a difficult  function  to  estimate 
from  the  spectra  themselves.  We  examined  a relaxation  modulation  on  a 
constant  ra^io  and  decided  that  T = 0.3  sec  was  an  appropriate 

relaxation  time.  Whatever  the  model,  the  correction  at  low  frequencies 
is  always  near  1.2. 

Finally  we  examined  the  receiver  crustal  response  and  assumed  an 
average  gain  of  2.3  for  the  vertical  component  of  the  P aves  and  2.5 
for  SlI  curves. 
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Then  at  the  long  periods  appri  prlate  for  estimating  r.ejsmlc  moment 
the  total  correction  is  of  the  order  of  the  geometric  attenuation  cor- 
rection (''<  5 X 10  ) , with  the  other  corrections  dividing  out.  At  high 
frequencies,  the  C)  correction  adds  a factor  o'f  up  to  20  (for  P waves) 
to  the  gain  which  must  be  applied  to  the  station  spectrum  to  get  a 
source  apectrum. 
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FIGURE  CAPTIONS 

Figure  4.1:  General  flow  diagram  for  autospectral  estimation  of 

seismogram  spectre. 

Figure  4.2:  Autocorrelation  lag  winaows.  Curve  1 is  the  rectangular 

window;  2 is  a Hanning  window;  3 ic  a Hamming  window; 

4 is  a Bartlett  window;  5 is  a Gaussian  window  with 

(1)  » 0.2;  and  6 is  an'jther  Gaussian  window  with  u - 0.J.11. 

Figure  4.3:  Average  autopowerspectral  levels  compared  to  the  perlodogram 

for  BAG  SPi:  (cop)  and  BAG  LP7  (bottom). 

Figure  4.4:  Window  closing  experlracut  DAG  SPZ.  TIk;  number  of  lags 

varies  from  5%  of  sequet  .-.e  ..ength  at  the  be  ttom  to  25Z 
at  the  top.  The  Fourier  anplltude  spectram  eatiaated  by 
direct  application  of  t.'if.  i'll  algortthm  in  shown  at  tlic 
top  for  comparison. 

Figure  4.5:  Window  closing  experiment  for  BAG  LPZ.  The  number  of 

lags  varies  from  102  of  sequences  length  at  the  bottom 
to  502  of  sequence  length  at  the  top.  The  Fourier 
amplitude  spectrum  estimated  by  direct  application  of 
an  FFT  algor ithia  is  uhowu  at  the  top  for  comparison. 

Figure  4.6:  Shapes  and  relative  amplitudes  of  the  four  waves  which 

contribute  to  the  telescismlc  P wave  from  a crustal 
source.  The  station  is  at  45°  distance,  40°  relative 
azimuth  from  a vertical  strike-slip  source  with  due 
North  slip  at  20  km  depth  in  a 30  km  thick  crust. 

Figure  4.7:  Total  P-wave  anelastic  attenuation.  The  highest  attenuation 

is  given  by  an  unmodulated  constant  ratio.  The 

modulated  attenuation  turves  arc  tangant  to  the  unmnciul.ited 
curve  at  the  rel..xation  frequency. 
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ABSTRACT 


If  the  specific  fracture  energy  associated-  with  an  earthquake 
rupture  is  non-zero,  then  the  rupture  velocity  cannot  exceed  the 
Rayleigh  wave  velocity,  , for  the  inplane  shear  node  of  extension 
or  the  shear  wave  velocity,  3 > for  the  antiplane  shear  inode  of 
extension.  The  radiation  efficiency,  > which  is  the  fraction 
of  the  available  energy  that  is  radiated,  depends  on  the  rupture 
hislrny.  The  radiation  efficiency  is  zero  if  the  crack  grows 
quasi-statically  and  improves  with  increasing  rupture  velocity.  The 
radiation  efficiency  approaches  unity  if  the  rupture  runs  nt  C for.  thr. 

Jv 

in)'1  .ne.  shear  i!.:ce  -i;u!  .iL  P for  lit;*  eut. ipla’Ki  shear  inode. 

If  a serai-infinite,  nntiplane  shear  crack  tip  accelerates  rapidly 

to  a terminal  velocity  v and  runs  nt  this  velocity  until  it  in 

arrested  abruptly  by  a fracture  energy  barrier,  then  tlie.  specific 

-2 

surface  energy  (ergs  cm  ) will  be, 
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where  L is  the  length  of  the  rupture  (cm) , is  the  trnctiori 

—2  ~2 
(dyne  cm  ) and  y in  the  rigidity  of  the  material  (dyne  cm  ) . Thus 

the  formula  for  the  specific  fracture  energy,  associated  with  the 

"seismic  gap"  model  of  fracture  arrest  (Hussei.nl  et  al.,  .197.')) 


Y 


o 


•• 


will  still  bo  valid  (to  v/lthin  a factor  of  tv.'o  for  v/fl  bcl; 'es 
•I* 

0 and  0.9  3 ) regardless  of  the  arrest  inc-.chanism. 
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